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ABSTRACT 

Paleomagnetic results from 5 stratigraphic sequences 
are reported. four of these concern the time interval 12,000 
to 13,000 years B.P. The motivation behind them was two- 
fold; first to seek evidence corroborating the existence of 
as extremely young geomagnetic polarity reversal - the so- 
called Gothenburg event, and second to explore the 
feasibility of using paleomagnetic signatures for strati- 
graphic correlation. No evidence for a geomagnetic reversal 
was found and evidence is growing that this event, if it 
existed at all, was more likely a non-dipole excursion of 
limited regional extent rather than a global polarity 
reversal. An enhanced clockwise secular variation loop was 
observed, and it is felt that this is indicative of some 
unusually intense non-dipole feature drifting westward past 
the sampling locality. This loop was identified at two 
localities some 250 km apart and the resulting correlation 
is of considerable significance geologically in that it 
establishes a stratigraphic relationship between two 
important volcanic tephras which had previously been 
unknown. 

The fifth section involves a detailed study of a single 
9,000 year paleomagnetic sequence. No geomagnetic excursions 
are found, and the absence of the Lake Mungo excursion 1s 
regarded as evidence supporting the spatially restricted 
nature of such phenomena. 


Spectral analysis (using the MEM technique) reveals the 
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presence of two types of periodic looping motion in the 
Riggins Road data. Highly elliptical counterclockwise 
looping with a period of roughly 5,000 years probably 
reflects a counterciockwise “wobble" of the dipole axis as 
suggested by other workers. The circular clockwise looping 
motion of period 1,900 years indicates the presence of 
westward drift of the non-dipole field some 30,000 years 
ago. 

The complex MEM technique is investigated by means of 
synthetic data sets in an effort to assess the validity of 
results obtained from very noisy records. Some interesting 
conclusions are obtained, the most important of which is 
that it enables one to identify both the frequency and sense 
of looping signals in the presence of noise levels twice 
that of the signal amplitude which makes it an ideal tool 
for analysing the type of time sequences encountered in 
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CHAPTER 1 
Introduction 


1.1 The Geomagnetic Field 


This thesis deals with the presentation, analysis, and 
possible implications of paleomagnetic data obtained from 
several young (< 40,000 yrs.) geological sequences. Such 
being the case, a brief overview of what is currently known 
(or hypothesized) about the earth's magnetic field and its 


behavior with time is first presented. 


To a very good approximation, the geomagnetic fieid can 
be represented by that of a magnetic dipole situated at the 
earth's center. The dipole axis presently intersects’ the 
Surface at 78./°9 'N, 289.5° E, about 11.5° from the rotation 
axis, and has remained in this position for as long as 
detailed world wide geomagnetic data have been available 
(roughly 150 years). Paleomagnetic data however shows’ that 
over much longer time periods (thousands to millions of 
years), the position of the geomagnetic pole has varied 
considerably. Motions taking place over millions of years or 
longer are generally called polar wander and play a very 
important role in global tectonic theory. Their time scale 
is much longer than the intervals studied in this thesis 
however. 

Changes in dipole orientation with periods of about 10%years 
also take place and these are often graphically referred to 


as "dipole wobble". There is as yet no widely accepted 
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theory as to the precise nature and period associated with 
this motion (Kawai and Hirooka, 1967; McDonald and 
Gunst,1968), but most authors agree that if averaged over a 
sufficient length of time the geomagnetic axis coincides 
closely with the rotation axis. In other words the time 
averaged geomagnetic field is thought to approximate that of 


a geocentric axial dipole. 


The actual degree to which the earth's field is dipolar 
is a matter of considerable importance in both paleomag- 
netism and geomagnetism since it allows one to specify a 
mathematical relationship between field direction at a given 
location on the earth's surface and the position of the cor- 
responding geomagnetic pole. This relationship is as follows 


(Irving, 1964); 


L' = arcsin(sin(L)cos(P) + cos(L)sin (P)cos(D) ) 
o* = > + arcsin(sin(P)sin(D)/cos(L*) ) 
(for cos(P) > sin(L)sin(L')) 
o* = © + 180 - arcsin(sin(P) sin(D) /cos(L*)) 
(for cos(P) < Sin(L) sin(L‘)) 
where Le = latitude of pole 
@*' = longitude of pole 
L = latitude of field location 
@ = longitude of field location 
D = magnetic declination 
and P = geomagnetic colatitude which is given by; 


2cot(P) = tan(T) 


where I = magnetic inclination 
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These formulae enable paleomagnetists to reconstruct 
the relative positions of the continents, and their 
associated lithospheric plates over geologic periods of 
time, and thereby provide much of the basis of the modern 
theory of plate tectonics. The same formulae are also used 
extensively by geomagnetists, both to study the motion of 
the pole itself, and as a means of comparing directional 
data obtained from widely separated locations. If sufficient 
time averaging has been achieved one speaks of "paleomag- 
netic poles", whereas if insufficient averaging has taken 
place one speaks of "virtual geomagnetic poles" (V¥GP*s) 


(McElhinny, 1973). 


Obviously if the earth's field were always perfectly 
dipolar VGP's obtained from geomagnetic directional data 
would vary only as a function of time and not as a function 
of the location from which the data were obtained. Analysis 
of world wide directional data obtained from the present 
field however, indicates that this is far from true. Figure 
1-1 shows that VGP's obtained from globally distributed 
sites typically vary 10° to 20° from the actual geomagnetic 
pole. The part of the field responsible for this deviation 


is called the non-dipole field. 


The dipole and non-dipole geomagnetic fields can be 
distinguished by spherical harmonic analysis of the entire 
field, the dipole field being represented by the lowest 


order terms in the expansion, and the non-dipole part being 
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represented by all higher order terms. The justification for 
regarding the fields separately is that they are thought to 
arise from two physically distinct processes in the earth's 
core. The dipole field is believed to be governed by a 
Single self exciting dynamo produced by means of large-scale 
fluid motion in the outer core (Jacobs, 1975). It is most 
likely driven by some form of convection and controlled to a 
large extent by those forces associated with the earth's 
rotation which would account for its axial nature. The non- 
dipole field, on the other hand, is thought to arise from 
small scale eddies or perturbations of the main current 
system close to the outer core boundary. Such features could 
be governed to some extent by either thermal Or 


topographical irregularities at the core-mantle interface. 


At the earth's surface the non-dipole strength is 
typically about 1/20 that of the dipole field (Irving, 
1964). The non-dipole field however is responsible for 
‘essentially all of the shorter period (tens to thousands of 
years) changes -- the so-called secular variation. Such 
variations have been observed directly for the last few 
centuries at various locations throughout the world. The 
magnitude of these geomagnetic directional fluctuations is 
typically of the order of 10° to 20°. fhe best known 
examples are those recorded at London and Paris (Fig. 1-2) 
which describe almost complete loops in magnetic direction 


with an apparent period of about 500 years. 
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The non-dipole field is generally displayed by means of 
isomagnetic maps on which one of the field elements (e. de 
declination, inclination, intensity) is contoured. The 
contours on such maps tend to close about regions of 
continental dimensions forming clearly visible cells several 
of which are present over the earth's surface (Fig. 1-3). 
Many of these features appear to grow and decay with 
lifetimes in the order of hundreds of years, whereas others 
appear to be stable over much longer periods of time. 
Perhaps the most interesting feature of the non-dipole field 
is its tendency to drift westward. Although observation of 
individual non-dipole features over the past century 
indicates that these features often follow guite different 
paths, and move at quite different rates with respect to one 
another, the overall tendency toward westward drift is 
unmistakable, the average drift rate being an easily obser- 


vable 0.189 per year (Bullard et al., 1950). 


Although few deny its existence, much concerning the 
actual nature of the westward drift remains unclear. The 
actual drift rate, for instance, is difficult to define, as 
the value obtained, even for a world wide average, depends 
on the method of analysis used. The value guoted above was 
obtained from calculations involving the non-dipole field 
elements directly. Analysis of the time derivatives of these 
same elements (the secular variation field) through 
comparison of isoporic charts for different epochs, gives a 


significantly higher drift rate (0.329 per year, Bullard et 
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ale, 1950). Yukutake and Tachinika (1968), attribute this 
difference to the presence of large stationary anomalies 
which tend to lower the mean drift rate calculated from the 
non-dipole fieid. Provided these anomalies do not fluctuate 
rapidly in intensity, their presence should not affect the 
velocity obtained through analysis of the secular variation 
field. Thus Yukutake and Tachinaka conclude that the higher 
the two velocities more truly reflects the velocity of the 
"drifting" portion of the non-dipole field. Simply dividing 
the non-dipole field into standing and drifting parts does 
not fully explain the marked regional variation in drift 
velocity. This is especially true in Canada and Alaska where 
very small, but definitely measurable, negative (eastward) 
drift rates have been observed (Cain and Hendricks, 1968; 


Yukutake, 1962). 


In addition to its present, rather complex, nature 
there is evidence to suggest that westward drift has changed 
with time as well. For example Harwood and Malin (1976) show 
that the rate of westward drift decreased by almost 50 
percent over the last few decades. Such a remarkably rapid 
change raises some doubt as to the long term stability of 
westward drift. Indeed, Hide and Malin (1970) have suggested 
that the drift may be a transient phenomenon initiated some 
500 years ago. In contrast to these observations, archaeo- 
Magnetic studies suggest that at least some features of 
Magnetic declination and inclination have probably been 


drifting westward for roughly 1000 years or more (Yukutake, 
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1962, 1967; Tarkhov 1965; Burlatskaya et al. 1965, 1968). 
Futhermore Bucha et al. (1970) have shown that a maximum in 
geomagnetic intensity may have been drifting westward for as 
long as 3000 years. Clearly there is considerable 
uncertainty regarding the long term nature of the secular 
Variation. Runcorn (1959), and Skiles (1970) have pointed 
out a meanS by which paleomagnetic data from a_ single 
location, can be used to infer the presence and direction of 
deyvit.. This technique, and its application to a particular 


set of paleomagnetic data is the subject of Chapter 3. 


Another geomagnetic phenomenon of considerable 
importance is that of polarity reversals. In contrast to the 
temporal Variations discussed above, these are dipole 
features and are associated with much longer time intervals. 
The last well established reversal occurred 0.69 X 106 years 
ago at the base of the present Brunhes polarity epoch. 
Recently, Morner et al. (1971) in Sweden claimed to have 
established a short reversal dated at approximately 12,000 
yrse BeP. and subsequently named the Gothenburg event. 
Corroboration of this event in widely separated geographic 
localities was necessary to establish it as a genuine 
reversal, and this was the original motivation behind the 
work described in Chapter 2. If verified, such extremely 
young and short events would be very important to geomag- 
netists in terms of the overall polarity spectrum and the 
mechanism of reversals, and to geologists as high resolution 


global time markers. Both these topics are discussed in 
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Chapter 2. It emerges that there probably are real geomag- 
netic events of this type but they do not represent genuine 
reversals. Currently they are given the label "geomagnetic 
excursions", several of which have been reported from around 
the giobe. They are thought to arise from sudden, rapid 
enhancements of particular non-dipole features in the 
outermost core and thus affect only limited geographic areas 
(Harrison and Ramirez, 1975). Verosub and Banerjee (1977) 
have recently provided a comprehensive review of the 
appropriate evidence for these features and conclude that 
the existence of any proposed excursion is not yet 
"sufficiently well established for its reality to be beyond 


question" (p. 153). 


1.2 Methods and Techniques 


The laboratory technigues involved in the measurement 
of the magnetic remanence vector corresponding to a given 
paleomagnetic specimen are now routine, and have been weli 
documented (see e.g. McElhinny 1973). MeasurementS were 
carried out using two commercial magnetometers, a Schonstedt 
model SM1 slow Spinner, and a Digico Balanced Fluxgate 
system (Molyneux, 1971). Specimens were partially 
demagnetized in alternating magnetic fields (A.F. 
demagnetization), using an apparatus built and described by 
Murthy (1969), in an attempt to remove unwanted percadecs 
conponents of magnetization. The actual fields used in 


demagnetization, are given later in the thesis as each set 
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of paleomagnetic data is discussed. 


Collecting of oriented specimens in the field presented 
some difficulty, as the material sampled was unconsolidated 
sediment. Two methods were used. One involved driving a 
brass channel horizontally into the vertical face of the 
outcrop; the other utilized small plastic cubes (available 
commercially) which were pressed into the sediment. Both 
techniques are described in detail in Appendix 1. A total of 
136 stratigraphic horizons from 5 separate localities were 
collected. The number of independently oriented samples was 
173, involving 225 individual specimens and a total of some 


1000 measurements of remanence. 


Results are analysed using standard Fisher (1953) 
Statistics, which can be regarded as the spherical analogue 
of Gaussian statistics. S. I. units are used throughout’ the 
thesis but since much geophysical literature is still 
written in cgs units, conversions for the most common 
parameters are given here: 


Magnetic Induction; 10-5 Gauss = 1 gamma 


1 nanotesla 


Intensity of Magnetization; 10-3 emu cm-3 = 1 amp/m 
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CHAPTER 2 
Late Quaternary Paleomagnetic Results from 


Southeastern Alberta and West Central Montana 


This chapter describes a paleomagnetic investigation of 
sedimentary deposits in the age range 12,000 to 13,000 
years, initially undertaken to seek evidence corroborating 
the existence of the so-called "Gothenburg Event", a 
supposed geomagnetic reversal recently reported from Sweden 
(Morner et al., 1971; Morner and Lanser, 1974). Recognition 
of this extremely young reversal would have been of 
considerable significance both geophysically and geo- 


Logically. 


A feasibility study was carried out on a section in 
southeastern Alberta. No reversal waS observed, but an 
apparently reliable geomagnetic record was obtained and it 
was decided to extend the study to other sections with two 
objectives in mind. First, to check the reality of certain 
distinctive features observed and thus improve our knowledge 
of long-term geomagnetic field behavior, and second to 
assess the possibility of using magnetic records as a_ means 
of regional geological correlation. Samples were therefore 


collected from three further stratigraphic sections. 


The four sections investigated are referred to by local 
geographical names -- Irvine, Choteau, Manyberries, and 


Diversion Lake (Fig. 2-1). Of these the first two are about 
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5 me thick and were sampled at more than 30 horizons’ each; 
the latter two are less than 2 m. thick, sampled at only 10 
horizons each. The Manyberries and Diversion Lake data are 
therefore of a reconnaissance nature, insufficent to stand 
alone but hopefully of some use in making comparisons with 
the more thoroughly investigated Irvine and Choteau sites. 
The vertical extent of each section nig indicated 
schematically in Figure 2-1 with the horizontal lines 
representing horizons sampled. The zero on the vertical 
scale represents the position of the Glacier Peak tephraat 
that locality. 


222 Geological and Sampling Details 


Three of the four sections studied contain a volcanic 
tephra originating from an eruption at Glacier Peak, 
Washington. This eruption has been dated at 12,750 + 350 
yrs. B.eP. (Mudge, 1967) using snail shells from the ash bed 
at Diversion Lake. The tephra was identified at all three 
locations on the basis of its chemistry (Lemke et. al., 
1975; Westgate, 1968), and is widely regarded as being the 
older of 2 ash-falls originating from Glacier Peak (Wilcox, 
1969; Westgate et. al.w,1970). However Westgate (Westgate and 
Evans, 1978) has discovered a third tephra in southeastern 
Alberta. This third tephra is undated and its stratigraphic 
relationship to the other two tephras is unknown, but it has 
a basically similar chemical composition and is therefore 
thought to have been erupted from Glacier Peak in the same 


general time interval as the other two ash falis. On the 
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basis of a chemical trend it seems likely that it is the 


oldest of the three (Westgate and Evans, 1978). 


At all four sections paleomagnetic samples were 
collected using a rectangular brass channel as described in 
Chapter 1. AS the geomagnetic features of interest in this 
study were relatively large compared to experimental error, 
it was decided to examine as many individual horizons as 
possible rather than concentrating on reducing noise by 
replicate sampling. To this end collection was restricted 
almost exclusively to one independently oriented sample per 
horizon. It was generally possible however, to cut and 
measure twWo or “more individual cubic specimens from each 
sample. Averaging the results from these specimens helps 
reduce noise arising from any distortion of the sediment 
during Sampling, and the measuring process in the 


laboratory. 


The Irvine section is Located in southeastern Alberta, 
approximately 10 km north of the town of Irvine (Lsd 4, sec 
25, fp 42, ® “2, .W 43 £50.19 Ny, 110239 WW). The sequence 
sampled is a near vertical, 5 metre thick, road cut exposure 
consisting entirely of horizontally bedded lacustrine silts 
(and some minor sands), probably of supraglacial origin. The 
marker horizon used was the undated Glacier Peak tephra 
discussed above -- at this locality it forms a_ white, 
clearly visible layer about 2 cm thick, located approximate- 


ly in the middle of the section. A total of 35 horizons were 


ov molgnne. ake 


se ing tame ab ry) ‘paanedo. ew 1 . +:  pnke ihe oe LL ran or 
ns AY gla ee ee BB hee | peo 


etl iv 1, 


pdt i teenasa to. dernsned of teat 


Meee 


Th ae 


mid 
ne op ae 


(aorre: tetnonheages of Soansn00 


os eitordod ‘oubivitne yom: 


ts onion padoabon ie caiertatonas eae yates ol tse0q 


bedpiagaes ; ‘a. sois2ehLes “t pie ae oe “wemhiqnea, gaso 1qe3 


wie} 


- req olyaae bet awixo + {susbangebmb one, poe ‘devie foe seonls 


Bix S00 ot svoniad BEERS ehtenaaoe Il 40 wos tied 


r =p 
oe3 aoat sae Loiges aide ieubawthat “stoi 40! owt. ‘eruesom 


eu) 
eo 


i i nme 


om tei. ssomioege inetd: ont avingert and) we 
“toonibiae. aay ho on sodas. xan, jue’ eatakes eaten 


eed. oe) nH 


edt) ones aeeneam Rokaienen, bh oan cao <a 
ay, Ad + he | i a vf ‘lh Nn, Seresmodal 
a oe a ay ull ean li ka 
ahh edits azote aed ven mt - hognoad eb aohtoee anas aa Ph 


‘ine dk er © ak 


Chel eoeg 


Mm 268 " best, oaivat to avo: de: Aeon ‘aa. ur ‘tdesonixe 
hi | f. pe a a dh 


Diccouiaed Sol oe e081 Pe 


us s! 3 i 


a H pePlinat old a 


18 


Sampled at 15 cm intervals, from 228 cm above the tephra to 


270 cm below it. 


The Choteau section is located some 250 km southwest of 
the Irvine section (NE 1/4, SW 1/4, sec 8, Tp 24N, R 3W, New 
Rockport Colony, Teton County; 47.89 N 112.0° W). The tephra 
at this locality is approximately 1 cm thick and lies 25 cm 
below the upper boundary of a thick sequence of glacio- 
lacustrine sediments. This sequence is overlain by a 
different depositional unit about 1.3 meters in thickness 
which contains the much younger Mazama tephra (approx. 6600 
ycse old (Lemke et ale, 1975)). The entire section is 
exposed on the nearly vertical wall of a 400 m. long 
landslide scarp. Lemke et al. describe this section in 
detail and indicate the presence of two depositional breaks 
in the lacustrine sediment below the Glacier Peak tephra. 
The poSitions of these depositional breaks are shown in 
Figure 2-2b. Lemke et al. argue that the lacustrine sediment 
beneath the lower depositional break was probably laid down 
during an earlier glaciation ( Bull Lake Glaciation), which 
implies an age greater than 70,000 years. The paleomagnetic 
results discussed below are of some importance in this 
context, Since they do not concur with this interpretation. 
Samples were collected at 40 horizons ranging from 11 cm 
above the Glacier Peak tephra to 468 cm below it. The 
average vertical distance between horizons was 15 cm, but a 
zone of interest was subsequently resampled on a more 


detailed basis. 
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The Manyberries section is located approximately 80 km 
southwest of the Irvine section, about 1 km south of the 
Village of Manyberries, Alberta (lsd 14, sec 13, Tp 5, R_ 6, 
W 4; 49,49 WN, 110.79 W). It consists of a meter thick 
sequence of horizontally bedded lacustrine sands, siits, and 
clays enclosing a horizontal band of white tephra approx- 
imately 3 cm thick. This sequence is part of a 10 meter 
vertical stream cut escarpment described by Westgate (1968). 
The lacustrine sediment overlies about 1 meter of sand and 
fine gravel with inclusions of till, and this in turn 
overlies a brown basal till. Above the lacustrine sequence 
lies a meter of late glacial sand and fine gravel, and 
alluvial sediment then extends to the top of the biuff. The 
tephra has been shown, on the basis of chemistry, to be of 
Glacier Peak origin, and identical to that found at 
Diversion Lake, and hence can be assigned an age of about 
12,700 years (Westgate and Evans, 1978). Paleomagnetic 
specimens were obtained from 10 horizons in the lacustrine 
seguence ranging from 60 cm above the tephra to 80 cm below. 
The vertical distance between horizons was typicaily about 


18 cm. 


The Diversion Lake section is located roughly 50 km 
southwest of Choteau, on the shore of Diversion Lake, Sun 
River, Montana (NW1/4, SW1/4, sec 36, Tp 22 N, R 9 Ws 47.69 
N, 112.79 W). Considerable digging into the side of the 


Shallow embankment at the edge of the lake exposed about 2 
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meters of lacustrine silts with interbedded sands and 
gravels and a band of white tephra about 2 cm thick. As 
discussed above, this tephra is thought to be the same as 
that Podaa near Manyberries:.and Choteau, and has an age of 
12,750 + 350 C-14 years (Mudge, 1967). Ten horizons were 
sampled from 26 cm below the tephra, to 154 cm above it. The 
presence of gravel beds in the sequence made uniform 
sampling impossible, however a typical distance between 
horizons was 20 cme A geological description of this 
locality is given by éalte et. al. (1975). 


223 Measurement 


and Results 


Remanence directions were measured using a Shonstedt 
model SSM-1 slow spinner magnetometer, Whenever possible (in 
roughly 80% of ail cases) 2 or more individual specimens per 
horizon were measured. Five pilot specimens were subjected 
to stepwise A..F. demagnetization in fields up to 180 mT. A 
typical median destructive field (MDF) was found to be about 
20 mit. In all but one case a relatively large angular shift 
took place on the first A..F. treatment (10 mT), suggesting 
the removal of a soft component of magnetization most likely 
picked up in the laboratory. Smaller shifts (roughly 10° ) 
were then observed for successively higher A..F. treatments, 
until eventually directional instability set in, anywhere 
between 30 and 80 mf. All remaining specimens were therefore 
demagnetized at 10 and 20 mf. A complete archive of ail the 
measurements iS given in Appendix 2. Table 2-1 summarizes 


the more important features. As with the pilot specimens, 


¥ Bt ate ar ao é 008 


eine “anoatved wt sleoer ; 
age ah ‘vedo c ser. hed vaniqa? 


P46 
any) see 


ekde | | ta” “ap haqisaaes | faoivol bey a Ce ae aay anowka ot ‘a 
> GaNREY of oo teed we novie a ‘whined Vayee 


ano3 x0 stem eonaupes: oh st ae 


saev ied, ‘eaaase li ‘Teotget 8 TeNen se 


: 
Sa af 7 
| hae samznemacaat ae Ee 
FY om Pe aay ui . 4 by on 7 
_Ahesanode 8 poten begat wtow “annkteersbb sosons tos, 70 a 


ak) af tdbexog. novensdW samagqnoteanm awaahgs sole penne abas 
or aneaigoqe benbe bone oten. $e a (avene “beens aa6 vidpudt |, 
hos see hata onew aitengooga, tebzq) ‘ake: er ood how kted a 

) ‘ait f oar ot qe shtek® a sosvectssapaaad iAsk subuqate ot : a | 
anode aa od: ‘banok aoe aan ogee 2M. ant & } aaebon a 
; tdde aekogns oumed eienaniag: —— sa ana te : 
 pedteonene tal OF a ates Pi i | 
F vienat seo, fot nxldo ages te | 
| ene videoony atte tie 


oh nu is 


(han i 


21 


average angular shifts between 10 and 20 mT were 
Significantly smaller than those between NRM and 10 mf for 
all four sections. For three of the sections within-site 
scatter eS iightre smaller at 10 mT than at 20 mT, but the 
Diversion Lake data is more highly grouped after a 20 at 
eataere: During laboratory storage (3 years) some of the 
Irvine specimens acquired a rather hard secondary 
magnetization which could not always be eliminated by 10 ni, 
In view of these observations, it was decided that the most 
appropriate treatments were 10 mf for the Choteau and Many- 
berries sections, and 20 mf for the Irvine and Diversion 
Lake sections. However it is stressed that this choice does 
not alter any essential feature of the conclusions drawn 
below. The results are illustrated in the form of dec- 
lination and inclination magnetograms in Figure 2-2 (a-d). 


2.4 Discussion 


In order to assess the statistical significance of the 
fluctuations in remanence direction apparent in the 
Magnetograms, the serial correlation test of Beran and 
Watson (1967) was applied to each of the four data 
sequences. AS Can be seen in Table 2-1, three of the 
seguences easily pass the test, indicating strongly that 
sequential geomagnetic field behavior has been recorded. The 
failure of the Manyberries sequence suggests that any geo- 
magnetic oscillations which may have occurred were of 
sufficiently small amplitude to be obscured by noise 


associated with sampling and/or the acquisition of magnetic 
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Paleomagnetic results from the 


a) 


b) 
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d) 


Lrvine section, Abberta 
(triangles = QV samples). 


Choteau section, Montana 
(asterisks = QP samples). 


Diversion Lake section, Montana. 


Manyberries section, Alberta. 


These magnetograms are plotted on different 


scales to provide maximum clarity compatible 


with page size. Stratigraphic elevations are 


positive (negative) above (below) each local 


tephra. Declination is: positive east of 


north, «and inelination, 15 positive couwnward: 
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Table 2-1 


Summary of Horizon Mean Data 


Irvine Choteau Diversion Manyberries 


Lake 

Stratigraphic 
thickness. (m) 4.98 4.449 Te REHe) LoAG 
Number of horizons 35 40 10 10 
Mean angular shift 

0-1LOmT IR NS: Sie oF hss LS Ray Ox 

10-20mT SENG 2 i ra EOwa 
Mean Qe 3 

NRM 6.6° 36) 2. Se 6702 

10mT 65291° ere dle. aia 5 oe 

20mt dea? Be OO MO MEA So Tee 
Mean Intensity 

NRM Graal TO aS) Ano 

'Cleaned' 2 i2 Gee ied 329 
Beran-Watson Test 

Serial Corre- 

lation Index Gao Lee Phas: URS 

Sig. Level >99% >99% >99% None 


Notes: 


-l 
Oe 3 is thevstandard error of the mean (- 81(KN) ~*) 
where k is Fisher's (1953) precision parameter and 
N is the number of specimens per horizon. 


2. Magnetic intensities are given in units of 107 2a/m. 


ry hes a sot We iuangiit! ato 


a | 
iy, aN 


Hy iat | 


Bait 


ary 


ae a 


Sm , a, 
sey : 4 a 
ee an i oe ae 
af aL 


eaten MO ee Ae ee oe gel aeteera . 
Pt es Ts Who on Ge at i aan 


7 io eo 
measoi aye 


28 


remanence. Nevertheiess acquisition of apparently meaningful 
geomagnetic records from 3 more or less contemporaneous 


sections made the prospects for correlation attractive. 


It was initially felt that the unusually large dec- 
lination swing discovered above the tephra at Irvine (Fig. 
2-2a) would be recognizable if found at other localities. A 
Similar feature is, in fact, visible on the Choteau 
magnetogram (Figure 2-2D). The comparison is better 
illustrated by plotting the appropriate data on a stereonet 
(Figure 2-3a). It is seen that an open loop having the same 
clockwise sense, similar magnitude, and occupying the same 
part of the stereonet is recorded at both sites. The reality 
of this loop is confirmed by the four QV horizons which 
represent a profile displaced about 3 meters laterally from 
the main profile and collected a year later. Some additional 
samples were also collected at Choteau (QP horizons) and 
they too are displayed in Figure 2-3a. The data points on 
each loop do not correlate on a one to one basis, however 
this is not surprising for two reasons. First, because of 
the distance between the two sites (250 km), slight 
differences in the secular variation patterns are to be 
expected, even in the absence of experimental error (e.g. 
see historical data from London and Paris, Fig. 1-2). 
Secondly, the sampling interval at Choteau is much finer 
‘than that at Irvine due tc subsequent resampling in the 
Vicinity of the loop. In particular, the five southerly and 


westerly points on the Choteau plot span less than 10 cm 
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Stereographic plots of paleomagnetic results 
from Choteau! €dous) anamirvine (CloriLamebes). 
Open circles correspond to triangles in 

Fie. 2-28 and asterisks: ai. bey. be 

a) Details of the secular vata ationwicoap 
(see text) recorded at both sites. 
phat data from =-394 vem ctor =2.06 vem, 
Levine) data rom +45 "em toy telioar cm. 

bo? Comparisonweor martes: §Ga) Susine. sa 
sliding 10 cm averaging window to smooth 
the densely sampled Choteau data. StraeiS 
graphic intervals ares Choteaus-330- cH 
toy-=Z250 cm: sivvine -105 vem toes) oem. 

c) Vectorial averages and their circles of 
95% contidence for the stratigraphic 
intervals immediately above and below the 
secular -variationwmloop, -)lustrated in. Cay) 
and (b). Stratigraphic intervals are civen 
Ln vem. The sequences of numbers in square 
brackets give the average declination and 
tnclination,. thesnumber of horizans 
involved, the precision parameter (k), and 


the circle of 957 scoptiidence in) despees. 
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whereas the sample spacing at Irvine is essentially 15 cm 
throughout. A more appropriate comparison between the loops 
is illustrated in Fig. 2-3b after the densely sampled 
portion of the Choteau data was smoothed using a sliding 10 


cm averaging window. 


Mean directional vectors for zones immediately above 
and below the declination swing at both sites are plotted in 
Fig. 2-3c. Magnetic vectors near the base of the Choteau 
sequence were very scattered (N=6, D=007, I=+55, k=8) but a 
subset of these produced a much better grouping (N=4, D=012, 
I=+49, k=36), and it is this latter group which is 
represented in the figure. Data at both sites appear to be 
considerably less eaaeecncd above the declination swing than 
below. Mean directions above the swing are not only statis- 
tically indistinguishable from one another, but also exhibit 
the same westerly trend (Figs. 2-2a, b), although lack of 


outcrop curtails the Irvine data. 


It is concluded that the loops recorded at the hee 
localities result from the same geomagnetic feature and are 
therefore contemporaneous, The geological implications of 
this conclusion are interesting. The fact that the loop at 
Irvine is above the tephra, whereas at Choteau it is below, 
proves that the tephras are not the same, and that Irvine is 
the older, as suggested by chemicai data (Westgate and 
Evans, 1978). Furthermore, unless one resorts to an 


inordinately large gap of some 50,000 to 60,000 years 
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between these two Glacier Peak ash falls, the correlation 
rules out the possibility that sediments below -240 cm at 
Choteau represent the much older Bulli Lake Glaciation as 


inferred by Lemke et al. (1975). 


The absence of any marked geomagnetic features at 
Diversion Lake (Fig. 2-2c), and the lack of any coherent 
sequential behavior at Manyberries limits the potential 
usefulness of these data. However as the chemical arguments 
for correlation of the tephra at Deyersion Lake, Many- 
berries, and Choteau are quite strong (Westgate and Evans, 
1978), it is expected that there should be some agreement 
between the magnetic results from each sequence. Compare 
EPCSE the data from Choteau and Diversion Lake. 
Unfortunately, geological conditions were such that nearly 
all samples collected at Choteau were from below the tephra 
whereas at Diversion Lake the majority came from above the 
tephra (Fig. 2-1). At Choteau, above the secular variation 
loop pee teeed earlier, directions are steadily west of 
north until the uppermost 30 cm or so are reached (Fig. 2- 
2b), where a rapid easterly swing takes place. Thus it is 
important to restrict comparison to the narrow zone of 
sampling overlap at the two localities. This is done in 
Figure 2-4. Despite the small number of horizons represented 
the overall agreement is reasonable, and if attention be 
restricted to the horizons immediately adjacent to the 


tephra the agreement is excellent. 
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Diversion Lake. Stereographic projection. 
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The Manyberries data is now considered. As at Choteau 
(below the tephra) declinations are mostly to the west and 
as the ash is approached an easterly swing is observed (Fig. 
2-2a), but dguectiens immediately adjacent to the ash do not 
agree with the Choteau - Diversion Lake results. In fact 
westerly decinaticns are reestablished above the tephra. It 
appears that either the Manyberries data has somehow 
suffered a systematic bias, or that the tephra at this’ site 
is not precisely equivalent in time to those found at 
Choteau and Diversion Lake, despite the chemical evidence. 
However, the Manyberries results are very scattered and fail 
the serial correlation test, therefore it is felt it is 


premature to draw any firm conclusions from then. 


No geomagnetic reversal is present during the time span 
represented by all four data sets. This is consistent with 
other recent evidence which casts considerable doubt on the 
reality of the so-called Gothenburg Event (Thompson, 1976; 
Sukroo et al., 1978) as a genuine geomagnetic reversal. It 
is also interesting to speculate on what sort of geomagnetic 
behavior gave rise to the feature observed at Irvine and 
Choteau. The fact that is takes the form of a clockwise open 
loop, iS reminiscent of the historical secular variation 
patterns recorded at many observatories. The loop discussed 
here however, is much larger. In the light of current 
knowledge concerning the geomagnetic field and secular 
variation it seems probable that this disturbance resulted 


from the rapid enhancement of some non-dipole feature in the 
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outermost core. The effect of such features can be of 
limited regional extent (Harrison and Ramirez, 1975) and if 
one were to drift west ward past a point on the earth's 
surface eee Chapter 3), a clockwise loop would be 
described, Similar to that observed in the paleomagnetic 


data discussed here, 


Finally, it is perhaps worthwhile to note that the 
above interpretation implies a time scale on the order of a 
few centuries for the observed geomagnetic loop to have 
occurred. This immediately puts constraints on the 
deposition rates at the Irvine and Choteau sections and such 
eo ey should be of considerable interest to geologists 
engaged in unravelling the Quaternary history of this 


region. 
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CHAPTER 3 
Late Quaternary Paleomagnetic Results 
from Southern British Columbia 


3.1 Introduction 


Since the beginning of the study of magnetization in 
recent sediments, interesting periodicities in the behavior 
of the geomagnetic field have been noticed (Johnson et al., 
1948; Nagata et al. 1949; Mackereth 1971; Creer 1974; Creer 
et al. 1976; Opdyke et al. 1972). It has, in fact, become 
generally accepted that the study of these periodicities is 
critical to the understanding of the physical mechanisms 
underlying geomagnetic secular variation. Variations in dec- 
lination and inclination have been examined in detail by 
many authors, but unfortunately in most cases the two 
directional elements have not been studied simultaneously. 
Where it has been possible to examine the two dimensional 
oscillatory behavior of paleomagnetic directions, a strong 
tendency towards looping has often been observed (Nagata _ et 
ale 1968; Denham 1974) as in the case of the well-known 
London and Paris historical records. The importance of this 
looping as a source of information concerning the drift of 
the non-dipole field with respect to the mantle, was first 
pointed out by Runcorn (1959) and is discussed in detail by 
Skiles (1970). In essentially all cases of practical 
interest the motion of the non-dipole focus past any point 
on the earth's surface will generate an elliptical loop in 


the geomagnetic vector at that point. In a normally 
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polarized field, the sense of the loop, looking along the 
vector, will be clockwise if the motion is westward, and 


counterclockwise if the motion is eastward. 


Previous studies have ne deea ay attempted to establish 
the presence of either looping or periodic variation in dec- 
lination or inciination simply from visual inspection. 
However, high noise levels, combined with the presence of 
more than one frequency in some data sets, renders visual 
analysis inadequate. In such cases, the application of more 
sophisticated methods of spectral analysis is desirable. 
Denham (1975), outlines the particular method of spectral 
analysis best suited to paleomagnetic data. The method 
invoives mapping the data onto the complex plane about the 
appropriate norm, and then analysing the resulting complex 
time series by means of the so-called maximum entropy 
method(MEM) of spectral analysis. The use of Burg's (1967, 
1968) maximum entropy spectral estimator was suggested 
because of its superiority over other techniques in dealing 
with short, relatively noisy data sets. The advantage of the 
complex representation of the data is that all the 
directional information iS contained in a single time 
series. As a result, one can determine not only the power at 
a given frequency, but also the tendency towards clockwise 
or counterclockwise looping at that frequency (for 


iilustrations see chapter 4). 


In this chapter, paleomagnetic data from a sedimentary 
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sequence in southern British Columbia are reported which 
constitute a time series of ancient geomagnetic directions 
Spanning some 9000 years. Fourier and the above MEM analysis 
is Ghaehednats and the implications of the resulting spectra 


are discussed. 


Paleomagnetic specimens were collected from a stream 
cut escarpment on the edge of Bessette Creek several hundred 
yards from Riggins Road approximately 10 km east of Lumby, 
Be Ce (50° i7* S5" N, 118° St' 45" WwW). The set of paleonag- 
netic results obtained from these specimens will henceforth 
be referred to as the Riggins Road data. The portion of the 
outcrop sampled (Fig. 3-1), consisted of dark, organic rich 
floodplain sediments containing four radiocarbon dated 
horizons and two tephras (Westgate and Fulton 1975). Two of 
the dates were obtained from detrital peat and the other two 
were from wood. The peat dates are not used in the 
calculation of deposition rate as both are considerably 
older than their stratigraphic position, relative to the 
wood dated horizons, would indicate. The tendency for 
detritus dates to be "too old" is not uncommon, whereas the 
wood, having retained its integrity, is less likely to have 
suffered a long pre-depositional history. Neither of the two 


tephras present have been independently dated. 


The lowermost wood dated horizon is at the base of the 
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section (31200 + 900 yrs B.P.). The younger horizon is 4.50 
meters above this and is dated at 25300 + 320 yrs BaP.«, 
giving an average deposition rate at 13 + 2 yrs/fcm. It is 
worth Reaewag out that the discrepancies between the 
various dates are not large -- a linear regression using all 
four dates yields a deposition rate of 10 yrs/fcm. Based on 
extrapolation of the estimated deposition rate, the age of 
the youngest specimens is 22100 yrs B.P. fhe oldest 
specimens were taken at the base of the section (31200 yrs 
B.Pe). With the exception of the two coarse sandy beds 
approximately 20 cm. in thickness near the top of the 
section, the sediment is remarkably uniform in character and 
there is no evidence of any major change in depositional 


environment during the time span sampled. 


Specimens were collected by pressing small plastic 
boxes into the vertical face of the outcrop. The boxes were 
Oriented using a Brunton compass-clinometer, then removed 
and tightly capped to prevent dessication. Two independently 
oriented specimens were collected from each of 37 4Strati- 
graphic horizons. The boxes were 2 cm. Square in cross 
section and thus each specimen spans some 2 to 3 decades; 
the horizon spacing is close to 19 cm which represents about 


250 years. 


The direction and intensity of the natural remanent 
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magnetization of each specimen was measured using a Digico 
balanced fluxgate magnetometer system. Eight pilot specimens 
were then subjected to stepwise A.F. demagnetization in 
fields up Le 40 mit peak value, and a typical median 
destructive field was found to be about 20 mT. The remanence 
directions varied by less than 10 degress until A.F. 
treatment exceeded 15 mT, beyond which some specimens became 
unstable. Consequently, 10 mT was chosen aS an appropriate 
field for A.F. cleaning of the remaining specimens. This 
treatment resuited in systematic movement away from the 
present field at the sampling locality, suggesting removal 
of small secondary components of recent origin. Prior to 
cleaning, intensities ranged from 2 x 10-! to 5 x 10-3 Aym 
with an average of 6 x 10-2 A/Ym. A complete list of results 
is given in Appendix 3. In figure 3-1 individual cleaned 
remanence directions, along with the vector mean for each 
horizon, are iliustrated in terms of declination and 


inclination magnetograms. 


The mean angular discrepancy between the pairs of 
specimens at each horizon is 5.99 + 4.49 (s.d.). Although 
this is the only measure of scatter that can be obtained 
directly from the data, a few more widely recognized and 
potentially more descriptive statistical parameters may be 
obtained indirectly, provided one makes some simplifying 
assumptions about the noise present. The statistical 
parameter generally used to describe the scatter in paleo- 


magnetic directions is Fisher's "k" (Fisher, 1953), given 
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Fige 3-1: Individual specimen directions and horizon means 
after A.F. cleaning at 10 mT. Ail significant stratigraphic 
features in the portion at the section sampled are 
indicated. The white regions, which make the up most of the 
stratigraphic diagram, at the left, consist of dark 
floodplain silt. 
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k= (Ne) Gi) 
where eivivs quakes of paleomagnetic directions, and R is 
their vector sum. Unfortunately this estimate is poor when N 
is small. Obviously one cannot use all 74% remanence 
directions directly to obtain k since directional variation 
between horizons is to be expected due to geomagnetic 
secular variation. An indirect estimate of k can, however, 
be obtained using information from ali 37 horizons if one 
assumes that the noise distribution at each horizon is the 
same. This Ne ots is reasonable when one considers the 
uniform character of the sediment being sampled throughout 
the section, as well as the fact that the same orientation 
procedure was used for all specimens. It is apparent fron 
the above furmula that to obtain Ky an estimate of R is 
required. To illustrate how this estimate is obtained, 
consider the following derivation where x(i), y(i)- and z(i} 
denote the 3 components of the ith remanence vector, V{(i), 


and V(i)eV{j) denotes the dot. product of the 2 vectors. 


5 N \2 /N a2 
R” =| )x(i); + pyc)" +) J2(i 
i Wh \i aA 


Pye) x (3) + cee + Tz (4) 24) 


ij ij ae) 
“8 | 

= Y) (x) 2G) + y(i)y(j) a z(i)2(5)) 
ij 


NN > 
Py v(i)ev(j) = N*<v(i) ev (5)> 
13 
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In the final line of the derivation <V(i)®V(j)> denotes the 
average of all possible dot products for all 74 remanence 
vectors assuming that all time dependent signals, which 
would Guucowneas variation between horizons, are absent. AS 
these signals will, in fact, always be present in the data, 
one cannot determine <V(i)eV(j)> exactly, however a reason- 
able estimate of this quantity can be obtained by 
calculating the dot product for each pair of vectors at each 
horizon, and then averaging the results for all 37 horizons. 
Having obtained this estimate, R is given by; 
R=Nsqrt (<V (i) eV (4) >) 
and therefore the best estimate of k can be obtained using 


the formula given earlier. 


Using the above technique, the best estimate of k for 
the data discussed here is 2441, The quantity analogous’ to 
the standard deviation for a normal distribution is given by 
theta63 = 81/sqrt(k) degrees, and is called the circular 
standard deviation (Irving,1964) 2. Given the above value for 
k, the circular standard deviation for the Riggins Road data 
is 5.229, which corresponds to a standard error of the mean 


of 3.69%. The parameters k and theta63 are, of course, only 


1 The standard two-tier analysis (Watson and Irving, 1957) 
yields k=246. 

2The circular standard deviation differs from the normal 
standard deviation in that 63% {as opposed to 68%) of all 
deviations from the "true value" are smaller than this 
quantity. 
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meaningful if the distribution of errors associated with the 
directional measurements 1s Fisherian (Fisher, 1953). 
Another statistical parameter often quoted in paleomagnetic 
literature, which is not based on any particular probability 
distribution, is the angular standard deviation given by 
arcos(R/N) (Wilson, 1959). For the present data, the angular 
standard deviation and the circular standard deviation are 
numerically the same to the accuracy quoted. The scatter in 
the data is most likely due to uncertainty involved in 
orienting the plastic cubes, combined with minor distortion 


of the sediment which may occur during sampling. 


The Beran and Watson (1967) statistical test for serial 
correlation (Epp et al. 1971) indicates that the probability 
of sequential ordering in the sequence of mean directions is 
greater than 99%, and this implies that at least some of the 
oscillations apparent in the magnetograms (E1Gs 3-1) 
represent real temporal fluctuations in the paleomagnetic 


field, the significance of which is discussed below. 


The overall mean direction for the entire sequence is 
D=005. 8°, T=+67.69 (N=37, k=329, alpha95=1.49 ) which 
differs significantly from the present geomagnetic vector at 
the eaapbenn Site (D=022.29, tI=+72.69 ). However, the 
observed mean inclination is virtually identical to that of 


an axial dipole (+67.5° ). The mean of the 37 virtual geo- 
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magnetic poles (VGP's) (Fig. 3-2) lies at 023.1° W, 86.3° N 
(N=37, K=125, A95=2. 1°, arcos (R/N) =7.2°), and is 
Statistically distinct from the earth's rotation axis. This 
pole, although "right handed", is not "far sided" which is 
of interest in connection with qidaones well known offset 
dipole model (Wilson, 1970, 1971, 1972). However, too much 
weight cannot be placed on this comparison because the 9000 
years represented at Riggins Road may be inadequate to fully 
average out secular variation and thus obtain ae reliable 
paleomagnetic pole (for a recent review of this problem see 


McElhinny and Merrill, 1975). 


The internal consistency, significant difference from 
the present local field, and high degree of serial 
correlation combine to strongly suggest that the Riggins 
Road data represent a reliable geomagnetic record. It is 
therefore interesting to note that there is no indication of 
any large angular shifts which could be correlated with the 
geomagnetic excursion reported from Lake Mungo, Australia 
(30,780-28,140 years B.P.; Barbetti and McElhinny 1972), of 
that from Mono Lake, California (25000-24000 yr B.P.~; Denham 
and Cox, 1971; Liddicoat and Coe, 1975). This either implies 
that some (or all) of the sedimentation and/or magnetization 
dates are in error, or that such excursions result from the 
spatially restricted effects of specific non-dipole sources 
as envisaged by Harrison and Ramirez (1975) and discussed in 
Chapter 1. Harrison and Ramirez show that if one models the 


source of such an excursion as a small magnetic dipole 
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ara geomagnetic poles (VGP's) for all 
3/7 or tZons. wlth theis mean cand ts circle 
of 95% confidence. LGRE inditeates, the 

pole caleulated from, che 965 International 
Geomagnetic Reference Field, and PF indi- 
cates the pole corresponding to the present 
field vector at the sampling site. The 
arrow on the perimeter indicates the 
Longitude of. the, sampling site. «Equal 
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located near the core-mantle boundary, then the regional 
extent of the magnetic disturbance at the earth's surface is 
surprisingly small (radius of influence is approximately 40° 
at latitude 309). Such an explanation is not consistent with 
the fact that Large apparently correlative excursions have 
been observed in the Gulf of Mexico (Freed and Healy, 1975) 
and in the Indian Ocean (Opdyke et. al., 1974). One 
plausible explanation which seems to be consistent with the 
facts, would be that these excursions are the result of a 
drifting non-dipole disturbance located along a line of 
latitude slightly south of the equator. Such a disturbance 
could guite possibly affect the three areas mentioned and 
yet have no effect at ilatitude 50° north (Riggins Road 


location). 


The absence of any unusual geomagnetic behavior, such 
as excursions, combined with a reasonably accurate knowledge 
of time spacing between specimens invites the use of some 
form of spectral analysis. Initially the declinations and 
inclinations were analysed seperately. Equal spacing between 
the horizons was assumed, and both the maximum entropy and 
fast Fourier transform techniques were applied. For both 
methods, the resulting declination spectrum is dominated by 
2 peaks, one at a period of about 2000 years and the other 
at about 900 years; and the inclination spectrun is 
dominated by a Single peak at approximately 5000 years. The 
maximum entropy and Fourier resuits are shown in figures 3-3 


and 3-4 respectively. In order to obtain an acceptable 
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spectrum using the Fourier transform technique it was 
necessary to apply a Hanning window to the data consequently 
reducing the effective length of an already short time 
series. pee wase of this drawback, and since MEM is regarded 
as a superior technique, especially where short, noisy data 
sets are concerned, it alone will be used in all subsequent 


analyses. 


The average spacing between horizons is 19.0 cme (250 
yrs), but 2 gaps greater than 25 cm. exist in the ecar a 
indicating that the assumption of uniform temporal spacing 
is only a fair approximation at best. In order to improve 
the initial spectral estimates, the data was therefore 
interpolated to a uniform sampling interval of 19.0 cm., 
prior to spectral analysis (interpolation was performed on 
each cartesian coordinate using natural cubic’ spline 
functions (Greville 1967)). The resulting spectra are also 
shown in Figs. 3-3 and 3-4. The "900 year" peak in the dec- 
lination spectrum appears much smaller after interpolation, 
leaving only one dominant peak at 2000 years. The single 
peak in the inclination spectrum has been somewhat enhanced 
by interpolation. The reader should bear in mind at this 
point that the plots in figures 3-3 and 3-4 are power 
density plots, hence power is equivalent to the area under a 


peak rather than its height (Lacoss 1971). 


The presence of a spectral peak at 2000 years is 


intriguing for two reasons. Firstly, oscillations of Similar 
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period have been previously noted by other authors (Creer et 
al 1976; Yukutake 1962); secondly, 2000 years is close _ to 
the time associated with one full cycle of the so-called 
westward amie of the non-dipole field using the mean drift 
rate calculated from historical data, (Yukutake 1962; 
Bullard et al. 1950; Adam et al. 1964). If this 2000 year 
periodicity is indeed a product of westward drift then it 
should be associated with clockwise looping of the geomag- 
netic vector as described above. Obviously some knowledge 
concerning the presence or absence of looping in the Riggins 
Road data is of fundamental importance, and the analysis is 


therefore pursued further. 


Denham (1975) has suggested a method of spectral 
analysis which yields information concerning the degree and 
sense of looping, in addition to the power, at any given 
frequency. Briefly the method involves analysing both the 
declination and inclination information simultaneously as a 
complex time series. Spectra obtained from complex series 
are often asymmetrical in that a peak in the "positive 
frequency" half of the spectrum is not the same size as its 
counterpart in the "negative frequency" half. The greater 
the degree of asymmetry at a given frequency, the more 
Circular is the looping at that frequency. In general, 
looping at arbitrary frequency, f, 1s elliptical with the 
length of the semi-major axis given by; 

| sqrt (P(f£))+sqrt(P(-f)) { 


and that of the semi-minor by; 
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| sqrt (P(f))-sqrt(P(-f)) | . 
where P(f) denotes the power associated with a spectral peak 
centered on f. The sense of the looping depends on the sign 
of P (£) -P (-£). A more complete description of the technique 


can be found in Denham's (1975) paper. 


In order to generate the appropriate complex series the 
sequence of interpolated directions was first mapped onto an 
angular equidistant polar projection centered on the 
vectorial mean of the data. Each direction was then 
converted to a complex number (x+tiy) using the transfor- 
mations x=rcos®, y=rsind, where r=distance {in degrees) fron 
the center of the projection 8=azimuth (measured clockwise 
on the polar projection), and i=sqrt(-1). As the maxinum 
entropy routine used to compute the spectra shown in Fig. 3- 
6 (Anderson 1974), was originally designed to analyse real 
data only, some minor modifications were necessary before 
analysing the complex time series. A FORTRAN listing of the 
routine used in analysing complex data for this chapter and 


chapter 4 is given in appendix 5. 


Before discussing the spectrum obtained from this 
analysis, it iS important to point out that the ability of 
the maximum entropy technique to provide a good spectral 
estimate for a given data set, relies on the choice of the 
correct number of terms for the prediction error filter 
(PEF). Analysis of the present data using many different PEF 


lengths showed that filters ranging from 11 to 20 terms gave 
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Similar spectra. Experimentation with synthetic data 
containing known frequencies (see Chapter4), but having a 
Similar length and signal/noise ratio as the present data, 
indicated Ae the best representation of the true spectrum 
was obtained using PEF iengths within the above range. A 
more objective criterion for choosing the correct filter 
length has recently been put forward by Ulrych and Bishop 
(1975), who suggest that the optimum filter length is that 
which minimizes what they refer to as the "final prediction 
error" (FPE).~ A plot of FPE versus filter length is given in 
figure 3-5. The position of the first minimum is somewhat 
surprising, since a two term filter would be much too short 
to resolve any appreciable detail. Application of this 
criterion synthetic data indicated that in order to obtain 
meaningful spectra the first minimum in the FPE had to _ be 
ignored, however the next most pronounced minimum did, in 
fact, seem to be a reliable indicator of the correct filter 
length. The second most pronounced minimum in Fig. 3-5 
appears at 13, hence this was the PEF length used in the 
final analysis of the data. The resulting spectrum is given 
in Fig. 3-6 and is clearly dominated by 3 major peaks. One 
of these strongly suggests clockwise looping at a period of 
1900 years, whereas the other two are indicative of 
elliptical counterclockwise motion at a period of about 5000 
years. In the shorter period region of the spectrum the 
power present is spread over a wide range of frequencies, 


and neither clockwise nor counterclockwise motion appears to 
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dominate. 


The spectral peaks at 5000 years appear reasonably 
Sharp, but their actual widths cover an extensive range of 
periods conservatively estimated as extending from 7500 
years to 4000 years. In other words, the period of 
oscillation represented by these peaks is not significantly 
different from the 6000 year period discovered by Opdyke et 
al (1971) in sediments from the Aegean Sea. The similarity 
becomes more convincing when one considers that in both data 
sets the oscillation is primarily a feature of the 
inclination. These authors suggest that this periodic 
variation may be due to a counterclockwise circular motion 
of the geomagnetic dipole axis. Our data definitely imply 
counterclockwise motion of the pole, but of about half the 


12° amplitude suggested by Opdyke et al. (1971). 


The spectral peak at 1900 years covers a range of 
periods from about 2100 years to about 1800 years. The high 
degree of clockwise power associated with this peak strongly 
Suggests that it represents the period associated with one 
full cycle of the westward drift of the non-dipole field. If 
this interpretation is correct, then it implies that at 
least some features of the non-dipole field are relatively 
stable for periods of 2000 years or more, and that the 
phenomenon of westward drift was in existence some 30,000 


years ago. 
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CHAPTER 4 
Spectral Analysis of Synthetic Data 


4.1 Introduction 


In the previous chapter MEM analysis revealed three 
spectral peaks in the time sequence from the Riggins Road 
section (Fig. 3-6). The amplitudes (obtained by numerical 
integration) and. frequencies (f) are 2.7° at f£=-0.05, 1.89 
at £=0.05, and 1.69 at £=0.14. Considering the rather high 
noise level (standard error of the mean = 3.6°) one is 
justified in questioning the ability of the maximum entropy 
technique to yieid reliable information concerning this 
motion. An experimental MEM study by Chen and Stegen (1974) 
using a real time sequence of length 24, containing a single 
sinusoid and random noise, revealed that "acceptable" 
spectra could not be obtained for noise levels greater than 
80%. In view of this rather discouraging result it seems 
that if one is to place any confidence in the "2000 year" 
spectral peak discussed earlier, it is essential to enquire 
whether or not the maximum entropy technique can 
consistentiy identify low amplitude motion in the presence 
of noise similar to that encountered in the Riggins Road 


data. 


In this chapter sequences of synthetic paleomagnetic 
data are generated containing superimposed looping motions 
of known freguencies in addition to random noise. Spectral 


analysis is then performed on each sequence in order to 
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determine under what conditions "acceptable" spectra can be 


obtained using the maximum entropy technigue. 


The synthetic data are generated by a FORTRAN IV 
program written by the author, called "wrtfreq". Details 
concerning user input are provided along with a listing in 
Appendix 4. Basically the program produces an elliptical 
periodic signal of a given frequency, simply by incrementing 
the following function: 

F(t) =Aexp (2qi ft) +Bexp (-27ift) 
where t=time, f=frequency, i=Sgqrt(-1), B=clockwise amp- 
litude, and A=counterclockwise amplitude. The semi-major and 
semi-minor axes of the resulting ellipse are given by A+B 
and jA-B{ respectively. The sense of the ellipse is  clock- 


wise if B>A and counterclockwise if B<A. 


Tie wordcis, tLovappLoxitliate sthe types of notse, appropriace 
to magnetic remanence vectors each vector in the synthetic 
seguence is rotated from itS original direction by a 
randomly chosen angle 9. Care is taken to ensure that 989 is 
chosen from a noise distribution which is, to a sufficient 
approximation, Fisherian. Assuming Fisher's "k" is not too 
small (k>3 is sufficient) then the probability that a given 
direction will be observed at an angle greater than ®' from 
the true direction is given by the following formnuia (Watson 
and Irving, 1957). 

Prob (€>0') =exp{-k(1-cos (6°) ) } 


A subroutine named GGEXP currently exists in the IMSL 


og 


we MATEO | a) ae ba ' : " | a 
wb raded aM es “f : 

oe potent: » ao art 0 
ae renter io 


aasiati i te eae roecapan sete * 20: fan 
HP ie “peasant paivetion ea 


5 


teem (asaya jasc a a fey 

ae ianalobess ‘ae hapa  genbergettst | eal” oh a 
Sut wollen nae ont Jobustigan an hedoetossnanen=k ban oom” 
av. tid. 6. oe weiitttn ithe one ‘Bie aaa soa. aAkave - 7 ; 
“tseite vad saaetay nae: fia ontion sebraiieel yaenybia : cup 

Ly en Be | | asoaer ow “tala aca Peay! ae i 


_prasonoanny exon ry ost? an opnwteerig 
Sahota ye ob ee saloon hoes ping 
mr aoarronty Apangiae. bid veda. bec 
ae a sa paanay. seaien 
~ natakanet & t, ee, aki. ty 
lhcags 70H. ae ‘ne mtaeaod i 
i ida « pa capaniaiiiiad on 


60 


library! which generates random deviates "x" with the 
following distribution: 
Prob (x>x')=exp (-x'/m) 

where m is a constant specified by the user. Using this 
routine the appropriate angular deviations are generated by 
setting xm equal to 1/k and choosing 9 from the formula 
§8=cos~!(1-x). After rotating a given vector away from its 
original direction by an amount 9 it is then rotated 
azimuthally by another randomly chosen angle 4. In order 
that the noise distribution be azimuthally symmetric 9 is 
chosen from a uniform noise distribution in the interval 


(0,27), using the IMSL routine GGUB. 


Long sequences (10000 points) of pure noise (no 
periodic signal present) have been generated using the 
procedure decribed above and submitted to the "chi squared" 
test for “a ‘Fisher distribution (Watson and “Irving, *1957). 
The results indicate that the noise is indeed Fisherian at a 
confidence level greater than 97%. It 1s important to 
mention however, that Short sequences (40 points) of noise 
produced by this routine sometimes deviate significantly 
from a Fisherian distribution. Fortunately the IMSL noise 
generating routines all require as input a "seed" integer 
(range 171 to 2147483647) which uniquely determines the 


"noise" output. In other words the subroutines will always 
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produce the same string of "random" numbers if the same seed 
integer is provided each time one of the routines is called. 
The above problem can therefore be effectively eliminated by 
using only pre-tested" seed integers; that is, seed 
integers which have previously resulted in the production of 
Fisherian noise. After generating and testing approximately 
50 noise sequences of length 40 using randomly selected seed 
integers, 10 of these integers were found to produce noise 
which is Fisherian at a confidence level of greater than 
75%. These 10 integers have been retained, and utilized in 
the production of all synthetic data sequences studied in 
this chapter (these integers are listed immediately 


following the listing cf wrtfregq in Appendix 4). 


fo iilustrate the type of periodic motion implied by 
various kinds of spectra, several relatively clean (low 
noise level) sequences of data containing looping 
frequencies of varying ellipticity, have been generated and 
analysed using the complex MEM technique. This exercise also 
serves as a useful test against possibile programming errors 
in the MEM routine, as the true spectra of such clean data 


are eaSily predictable. 


The resuiting power density plots are presented along 
with polar projections of the 4 synthetic data sets in 


figure 4-1. The complete absence of counterclockwise power 
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Fig. 4-1: MEM spectra for signals of varying ellipticity 
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in the spectrum associated with circular clockwise looping 
is demonstrated in figure 4-la. The spectrum of a 2:1 clock- 
wise ellipse is shown in figure 4-1b where clockwise power 
outweighs counterclockwise Povwemy iby ana ehactor, i6fius9.. The 
spectrum in figure 4-1c represents that of a 3:1 counter- 
clockwise ellipse (ratio of counterclockwise to clockwise 
power is 25 to 9). Finally, figure 4-1d shows the spectrun 
resulting from a simple iinear oscillation where, of course, 


counterclockwise and clockwise power are equal. 


en SS ee eee ee ee >So ee SS a ee Se 


As the noise in a paleomagnetic time sequence wiil 
often be large compared to the relevant signals present, it 
is important to have some indication of the limits under 
which the complex MEM technigue can perform successfully. In 
this section an attempt is made to determine the smaliest 
Signal detectable in a paleomagnetic time sequence of fixed 
length (n=40)1 containing a noise level typical of real 
data. The synthetic data sets originally generated for this 
purpose each contain 3 superimposed circular looping motions 
of different frequencies but of the same amplitude (angular 
radius). Two of the signals are of clockwise sense at 


frequencies of 0.05 and 0.31 and the third is counterclock- 


1 This is similar to the length of the Riggins Road data 
(n=37) . 
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wise at a frequency of 0.14. The presence of 3 frequencies 
as opposed to one in the synthetic data provides a check as 
to whether or not the "detectability" of a signal isa 
function of its frequency, as well as its amplitude. The 
angular standard deviation of the noise added to the data is 
fixed at (5/sqrt(2))°. This noise level is considered 
appropriate since 59 is a typical standard error for a 
Single paleomagnetic specimen, and 2 specimens per horizon 
is regarded aS an optimum sampling scheme (Watson and 
Irving, 1957). This noise level is also very Similar to that 
found in the Riggins Road sequence, where 2 specimens per 


horizon were collected. 


Ten distinct synthetic data sets were generated 
containing signals of amplitude 1.0° and noise of the level 
indicated above, using ten pre-tested seed integers. The 
Same ten seed integers were then used in the production of 
ten more data sets containing signals of amplitude 1.59. 
This procedure was continued until a total of 50 data sets 
were produced with signal amplitudes of 1.09, 125°, 2.0°, 
2.59, and 3.0%. The power density spectra of the first 5 
data sets at amplitudes 1.09%, 2.0°, and 3.09 are shown in 
figure 4-2. One can see that all of the spectra shown at 
amplitudes 2° and 3° are reasonably good whereas 4 out of 5 
of the spectra at amplitude 1° are definitely unacceptable. 
It should be mentioned at this point, that in obtaining ail 
of the spectra discussed in this section, and henceforth in 


this chapter, the FPE criterion was utilized in choosing the 
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Fig. h—25 Power density spectra for synthetic data 
containing signals of varying amplitude and a constant noise 
‘level (see text). NR"n" denotes the particular noise 
realization at a given amplitude. 
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optimum prediction error filter length. It is important to 
note however, that in using this criterion on a large number 
of synthetic data sets similar to those described above, it 
was found that occasionally an FPE minimum would occur at a 
filter length which was far too low to produce a spectrum 
with acceptable resolution. An acceptable spectrum could 
generally be obtained in such a case by using the filter 
length where the next FPE minimum occured. Since PEF lengths 
of less than 6 were found never to produce spectra of 
sufficient resolution when analysing noisy data, an 
arbitrary lower limit of 6 was imposed each time the FPE 
criterion was used. An upper cutoff of n/2 (n=number of data 
points) was also imposed as suggested by Ulrych and Bishop 
(1975) « The filter lengths (FL) used in calculating spectra 


for the above synthetic data sets are given in Table 4-1. 


An obvious method of judging the “acceptability" of 
synthetic spectra, such as those in figure 4-2, is_ to 
compare the heights of the "true" spectral peaks with those 
of the ever present "spurious" peaks!., The power density 
ratio of the smallest true peak to the largest spurious peak 
is therefore chosen as a numerical indicator of the quality 


of a particular spectrum. Such a ratio will be referred to 


1 If a prediction error filter length of M is used when 
employing the MEM technigue, M-1 peaks will be present in 
the resulting spectrum, regardless of the true spectral 
content of the data. 
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here as an RTS (ratio of true to spurious) value. RTS values 
for all 50 synthetic data sets are represented graphically 
in figure 4-3). These values are also listed in Table 4-1. 
When examining this data one can generally regard spectra 
with an RTS value of less than 2 to be unacceptable since an 
observer who is not aware of the number of "true" peaks 
present (as is the case with real data) will find such 
spectra very difficult to interpret. This being the case, 
figure 4-3 indicates that attempting to identify signals of 
amplitude less than 2° using the complex MEM technique in 
noise of the level indicated above, would not be advisable. 
One should bear in mind however, that any conclusions 
concerning the limitations of the complex MEM technique 
drawn on the basis of the above results, may be overly 
pessimistic due to the requirement that 3 signals be 


identified at a given amplitude rather than one. 


As mentioned earlier, the presence of 3 signals in the 
synthetic data provides a rough check as to whether of not 
the rate of signal deterioration with decreasing amplitude 
is influenced by frequency. A slight trend towards improved 
Signal recovery at lower frequencies is, in fact, apparent 
but the effect of this trend on the observed power density 
is small compared with the effects of varying filter length 


and signal to noise ratio. 


In an attempt to assess the performance of the complex 


MEM technigue under less adverse, and perhaps more realistic 
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Fig. 4-3: Semi-log plot of RTS versus amplitude for 
synthetic data containing 3 superimposed periodic signals. 
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Table 4-1 
Summary of Results for Synthetic 
Data Containing 3 Periodic Signals 
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conditions, another series of synthetic data sets were 
. generated identical to the first, except this time 
containing only one periodic signal, at a frequency of 0.14. 
The results are Summarized in the same manner aS _ above; 
first graphically in figure 4-4, and then in table form 
(Table 4-2). Examination of figure 4-4 reveals that in 
typical paleomagnetic noise of the level described earlier, 
one can be reasonably certain of recovering a_ single 
periodic Signal with an amplitude as low as 1.59 In other 
words, the lowest allowable signal to noise ratio that can 
be present if one wishes to obtain a reasonable spectrur, 
given a data sequence of about 40 points, is approximately 


0.42. 


These results appear to conflict with those of Chen and 
Stegen (1974) who were able to obtain "acceptable" spectra 
in only 6 out of 10 runs with a signal to noise ratio of 
2252 The disagreement iS more apparent than real however, 
for 2 reasons. The first is that the time sequences used by 
Chen and Stegen were only 24 points long as opposed to the 
40 point sequences used here. The second, and probably more 
important reason however, is that the definition of 
"acceptable spectra" used by Chen and Stegen, although not 


clearly defined in their paper, appears to be considerably 
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Fig. 4-4: Semi-log plot of RTS versus amplitude for 
synthetic data containing only one periodic signal. 
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Table 4-2 
Summary of Results for Synthetic 
Data Containing One Periodic Signal. 
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more conservative than that used in this chapter?. 


The best estimate for the angular standard deviation 
associated with the noise in the Riggins Road data was found 
to be 3.6° Taking the amplitude of the "2000 year" 
oscillation discussed in Chapter 3 to be 1.6° gives a signal 
to noise ratio of 0.44. The analysis of synthetic data 
described above therefore implies that the "2000 year" 
oscillation observed at Riggins Road, is real but lies very 


near the lower limit of the complex MEM technique. 


One important conclusion drawn in Chapter 3 is that the 
power density spectrum of the Riggins Road data is 
indicative of counterclockwise elliptical motion of the geo- 
Magnetic vector at a period of approximately 5000 years. The 
assertion that the motion is counterclockwise rests on the 
fact that the counterclockwise amplitude associated with the 
Signal is larger than the clockwise amplitude by 
approximately 0.8°. Considering the noise level present in 
the Riggins Road data one is justified in challenging the 


significance of this difference. 
To obtain some indication of the fluctuations in 


i Spectra with RTS values as low as 2 have been generally 
regarded as “acceptable"., 
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apparent amplitude that one might expect in the spectra of 
noisy data such as that encountered at Riggins Road, 10 
synthetic data sequences were generated, each containing a 
single linear oscillation in addition to noise of the same 
level used in the previous section. In each seguence the 
total amplitude of the oscillation was 5° (ise. clockwise 
amplitude=2.59 and counterclockwise amplitude=2.5°) and the 
frequency was 0.05. Spectra were obtained for each of the 
data sets and numerical gieegnetien was performed to obtain 
the apparent clockwise and counterclockwise amplitudes 
associated with the ened present. The results are 


presented in Table 4-3. 


Considerable fluctuation is evident in the amplitudes 
shown, the average standard deviation being approximately 
15% of the true amplitude. Also evident is a slight tendency 
for the amplitudes to be too large. This trend is nost 
jikely a consequence of setting the integration limits too 
far apart when calculating the area under the spectral peaks 
(the limits used were .03 and .08 for both. positive and 
negative frequencies). The RMS deviation of the differences 
(clockwise amp. minus counterclockwise amp.) from their true 
mean of zero is 0.379, and the maximum value is 0.579%. As 
the difference value associated with the "5000 year" 
spectral peaks in the Riggins Road data is 0.89, it is 
highly unlikely that the apparent counterclockwise sense 
associated with this periodic motion, is purely a result of 


noise in the data. 
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Table 4-3 
Summary of Results for Synthetic 
Data Containing a Single Linear 
Oscillation of Amplitude 5.0° 
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As stated in section 4.3, given the noise level present 
in the Riggins Road time sequence, along with its limited 
length of 37 points, a Signal with amplitude equal to that 
of the so called "2000 year" oscillation found in the 
Riggins Road data, is probably the smallest Signal that one 
could hope to identify with reasonable certainty using the 
complex MEM technigue. One must keep in mind that section 
ue3 is designed to explore the limits under which the MEM 
technique can perform successfully, rather than to 
accurately simulate any real data set. For this reason there 
are several differences between the synthetic time sequences 
used in that section and the Riggins Road sequence. In view 
of the rather borderline nature of the "2000 year" spectral 
peak, these differences may be important enough to alter any 
interpretation as to the significancerof£ that peak, based on 
the resuits of section 4.3. In this section an attempt is 
made to produce synthetic data sets which simulate the 


Riggins Road data as closely as possible. 


Ten synthetic data sets of length 37 were generated 
containing signals similar to those thought to exist in the 
Riggins Road sequence (amp. 2.679 at freg.=-0.05, amp. 1.849 
at £=0.05, and amps 1.63° at £=0.14). The angular standard 
deviation of the noise added to the data was set at 3.6°. 
The synthetic sequences were also produced with the same 


unequal spacing in time (elevation) as the Riggins Road 
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data, and interpolation was performed prior to spectral 


analysis aS in the case of the real data. 


Spectral density plots for the ten data sequences are 
Shown in figure 4-5. Numerical integration was carried out 
to determine the amplitudes associated with each of the 3 
peaks, in each spectrun. These results, along with the RTS 
values and PEF lengths for each Sey are presented in 
table 4-4. If one judged the spectra purely on the basis of 
RTS value, then 8 out of 10 spectra could qualify as 
acceptable (RTS > 2). Most readers however, would likely 
agree on the basis of visual examination of the power 
density plots, that data sets 5 through 8 give unacceptable 
spectra. From these results it is concluded that, given the 
Signal to noise conditions present in the Riggins Road data, 
the chances of extracting, in recognizable form, the signals 
discussed in Chapter 3 through spectral analysis is between 


60% and 80%, assuming the signals are in fact real. 


One method of assessing whether the spectrai peaks 
found eh the Riggins Road data are real or spurious is to 
compare the spectral densities of the largest spurious peaks 
found in the synthetic data sets, with those of the "true" 
spectral peaks found in the real data. To carry out this 
exercise properly, it was necessary to recalculate the power 
density spectra for each data set using a PEF length equal 
to that used in analysing the Riggins Road data. This is 


because the power density of any spectral peak obtained 
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Power density spectra for simulated Riggins Road 
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Table 4-4 
Summary of Results for Simulated Riggins Road Data 
(unit of amplitude is degrees) 
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using MEM varies strongly as a function of filter length. 
Upon examination of the spectra for all 10 synthetic data 
sets calculated using a filter length of 13, only 3 spurious 
peaks (out of a total 120) were found to have power 
densities greater than 100 deg? per unit frequency, the 
largest of these having a height of 171 deg@ per unit 
frequency. The distribution of power densities for the 
spurious peaks is illustrated in figure 4-6. The shortest of 
the 3 peaks found in the Riggins Road sequence (the clock- 
wise peak at period of 5000 yrs.) has a power density of 269 
deg@ per unit frequency. LL one assumes that the 
distribution in figure 4-6 is tlog-normal!, then these 
results indicate that the probability of such a= peak 


occurring by accident is less than 1%. 


1 The power density values pass the chi-squared test for 
such a distribution (see appendix 6). 
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CHAPTER 5 
Summary and Conclusions 

Paleomagnetic results from 5 stratigraphic sequences 
have been reported. Four of these are discussed in Chapter2 
and concern the time interval 12000 to 13000 years B.P. The 
motivation behind this work was Penn eadia- first to seek 
evidence corroborating the existence of the so-called 
Gothenburg event (Morner et al., 1970), and second to 
explore the feasibility of using paleomagnetic signatures 
for stratigraphic correlation. No evidence for a geomagnetic 
reversal was found but this was not totally unexpected as 
more recent evidence (Thompson, 1976; Sukroo et alsa, 1978) 
indicates that the Gothenburg event, if it existed at all, 
was more likely a non-dipole excursion of limited regional 
extent rather than a global polarity reversal. An enhanced 
clockwise loop in remanence direction was however observed, 
and it is felt that this is indicative of some unusually 
intense non-dipole feature drifting westward past the 
sampling locality. It is tempting to equate this disturbance 
with that which initiated the Gothenburg "event" in Sweden, 
but until further evidence is presented supporting the 
existence of the Swedish excursion, such a conclusion could 
only be regarded as speculation. The existence of the loop 
was established at two localities some 250 km apart. This 
correlation was of considerable significance geologically in 
that it established a stratigraphic relationship between two 


volcanic tephras which had previously been unknown (Westgate 
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and Evans, 1978). 


In Chapter 3 a single 9000 year PMisaungnetes sequence 
is discussed in detail. No geomagnetic excursions are found, 
although two have been documented for the time span 
represented by the data (22000 to 31000 yrs B.P.). The 
absence of the Lake Mungo excursion at the sampling locality 
is regarded as evidence supporting the spacially restricted 
nature of such phenomena (Harrison and Ramirez, 1975). The 
absence of the Mono Lake excursion is somewhat more puzzling 
as it was discovered only about 1200 km from Riggins Road. 
It therefore appears that either a large depositional break 
is present in the section, or the dates at Riggins Road 
and/or Mono Lake are in error, or that the Mono Lake 
excursion did not occur at the time claimed. The latter 
hypothesis is supported by Verosub (1976) who was unable to 
detect the excursion in a contemporaneous seguence only 300 


km from Mono Lake. 


Spectral analysis (using the MEM technique) reveals the 
presence of two types of periodic looping motion in the 
Riggins Road data. The highly elliptical counterclockwise 
looping with a period of roughly 5000 years may be 
indicative of eastward drift of some portion of the non- 
dipole field, or of a counterclockwise -circular motion of 
the geomagnetic dipole axis as suggested by Opdyke et al. 
(1971). The circular clockwise looping motion of period 1900 


years indicates the presence of westward drift of the non- 
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dipole field. As mentioned in Chapter1l, there is presentiy 
considerable uncertainty regarding the long term stability 
of the westward drift thus any scien ee concerning its 
existence in past is of considerable importance. If one 
accepts the interpretation that the 1900 year period 
represents one full cycle of the drift, then the results of 
Chapter 3 not only demonstrate the existence of westward 
drift some 30000 years ago, but give an estimate of the 


drift rate for that time interval as well (0.19°/yr). 


In Chapter 4 the complex MEM technique was investigated 
by means of synthetic data sets in an effort to assess the 
validity of the results of Chapter 3. Some interesting 
conclusions were arrived at. The most important of these, is 
that the three principal peaks in the spectrum are not 
Simply products of the noise present in the data, and there- 
fore almost certainly represent real geomagnetic phenomena. 
Furthermore the sense of the periodic signals as inferred 
from the Riggins Road spectum is statistically significant. 
In spite of the above however, the findings of this chapter 
strongly suggest that the Signal to noise conditions under 
which the Riggins Road resuits were obtained are marginal 
(especially where the "1900 year" peak is concerned). In 
other words, anyone wishing to verify the presence of the 
above signals using a different paleomagnetic data _ set, 
would be well advised to either reduce noise by collecting 
more samples per horizeon, or collect a longer data sequence. 


Due to the marginal character of the Riggins Road Sequence, 
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the presence of shorter period (less than 1900 yrs.) 
oscillatory motion of the geomagnetic vector with an 
amplitude similar to, or less than that of the 1900 year 


oscillation, cannot be ruled out. 


Apart from the results concerning the Riggins Road 
data, this chapter has also demonstrated the impressive 
power of the complex MEM technique. The fact that it enables 
one to identify both the frequency and sense of looping 
Signals in the presence of noise levels twice that of the 
Signal amplitude, makes it an ideal tool for analysing the 


type of time sequences encountered in paleomagnetisn. 
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APPENDIX 1 
Description of Sampling Techniques 

Essentially all specimens from the four sites discussed 
in Chapter 2 were collected by driving a rectangular 3-sided 
brass channel into the vertical face of the outcrop. Care 
was taken to ensure that a sufficient portion of the chanrel 
was exposed after emplacement so that its orientation could 
be measured. Orientation was defined by three angles -- one 
called the "strike" was the horizontal bearing of the sides 
of the of the channel -- another, the "dip", was the angle 
of the top edge below the horizontal (measured in the 
vertical plane) -- the third, calied the "roll", was roughly 
equivalent to the rotation of the channel about its long 
axis, but more precisely was the dip of the top of the rear 
face (also measured in the vertical plane). A clockwise roll 
(iee. top of the rear face dipping to the right) was defined 
as positive. Measurements were taken with a Brunton compass- 


clinometer and were accurate to roughly one degree. 


Prior to removal of the channel from the sediment, a 
flat steel chisel was driven into the outcrop sliding along 
the top (open) face so that sediment in the channel would be 
enclosed on all four sides. Both the chisel and the brass 
channel were then gripped together and pulled from the 
outcrop simultaneously. To force the sample out of the 
channel a rectangular brass piston was placed behind the 
sediment and attached to a rod inserted through a hole in 


the rear face of the channel. Employing a device designed by 
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John Westgate, the channel was held in a stationary position 
while sufficient pressure was applied to the piston rod to 
force the sediment out. Once removed, the rectangular sample 
was marked to preserve its orientation then immediately 
sprayed with laguer to prevent it from drying out and 
crumbling. The samples were generally long enough so that 2 
to 4 individual cubic specimens could be cut from each. 
These specimens were sequentially labelied A, By, Cy, aa-« 
where A would be the one closest to the end of the sample at 


the outer face of the outcrop. 


As mentioned earlier the above method was used 
extensively in collecting the data analysed in Chapter 2, 
and although reasonably sucessful, the technique had two 
major disadvantages. Firstly the process of cutting the 
rectangular samples into cubes was extremely time consuming 
and sometimes impossible as the sediment would generally 
crumble if not handled properly. Secondly most samples were 
Significantly compressed while being extruded from the brass 
channel, causing distortion of the sediment sufficient to 
affect the magnetic remanence. Visual examination of samples 
in the field both before and after extrusion revealed that 
most of the compression was taken up in the centimeter of 
sediment closest to the brass piston. For this reason, 
remanence directions obtained from the so-called A specimens 
(see above) are not considered in the analysis of paleomag- 


netic data in Chapter 2. 
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In an attempt to avoid the above two problems a some- 
what simpler sampling scheme was devised to collect the data 
presented in Chapter 3. In place of the brass channel 
transparent plastic cubes, roughly two centimeters on a 
Side, were either pressed or gently tapped into the vertical 
face of the outcrop. Once in place, the vertical position of 
the sample was measured and recorded as in the previous 
method; the orientation procedure however, was somewhat 
different. The strike in this case was defined as the 
horizontal bearing of the rear face of the cube. The dip, 
although defined the same as previously, was measured by 
holding the Brunton flush against the rear face ina 
vertical plane perpendicular to the face. The roll was 
obtained in the same manner as before. Measurements were 
accurate to within 2° or 39°. After removal from the outcrop, 
excess sediment was carefully scraped away so that the 
sample was flush with the open side of the cube. The cube 
was then capped (each had its own tightly fitting lid), 
marked to preserve its orientation, and labelled with a felt 
pen. In the laboratory, the samples were simply measured 
while still in the cubes (the plastic cubes were totally 
non-magnetic) avoiding the weeks of work associated with 
cutting up the “brass channel specimens". As it was not 
necessary to force the sediment out of its container as was 
done previously, distortion wWaS Minimized as well. The one 
disadvantage of the newer scheme was that each oriented 


sample yielded only a single specimen for magnetic 
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measurement. For this reason at least two independently 
oriented samples were collected per horizon as opposed to 


only one when the brass channel was used. 
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The data used in Chapter 2 is listed here in the 


column 


colunn 


column 


column 


column 


column 


column 


QZIO1A 
QZIO1A 
QZIO1A 
QZ1I0 1B 
QZI01B 
QZ1I01B 
QZHO1A 
QZHO1A 
QZHOTA 
QZHO1B 
QZHO1B 
QZHO1B 
QZHO1B 
QYGO1A 
QYGO1A 
QYGO1A 
QYGO1B 
QYGO 1B 
QYGO1B 
QZGO1A 
QZGO1A 
QZGO1A 
QZGO1B 
QZGO1B 
QZGO01B 
QYFO1A 
QYFO1A 
QYFO1A 
QYFO1B 
QYFO1B 


1-- 


2-- 


12. 66 
12. 66 
12. 66 
11. 60 
11.60 
11. 60 
12.311 
12.31 
12.31 
11.68 
11.68 
11. 68 
11.68 
13.44 
13.44 
13.44 
12.31 
12. 31 
12.31 
12.32 
12.32 
12632 
12.06 
12. 06 
12.06 
13.55 
13:35 
ye re 
12.24 
12.24 


APPENDIX 2 


following seven-column format: 


specimen name 


specimen mass 


vertical distance above tephra in cm 


A..F. 


in grams 


magnetic moment in amp m2 


declination 


inclination 


Irvine Section 


228.9 
228.0 
228.0 
228. 0 
228.0 
228.0 
210.0 
210.0 
210.0 
210.0 
210.0 
210.0 
210.0 
195.0 
195.0 
195.0 
195.0 
195.0 
195.0 
180.0 
180.0 
180.0 
180.0 
180.0 
180.0 
165.0 
165.0 
165.0 
165.0 
165. 0 


1. 445E-07 
6. 460E-08 
3.446E-08 
8.961E-08 
5. 431E-08 
32871E-08 
1. 393E-07 
32597E-08 
22-508E-08 
72 T46E-08 
3. 696E-08 
2e 728E-08 
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QZM0O1B 
QZM0O1B 
QZMO1C 
QZMO1C 
QZMO1C 
QYMO1A 
QYMO1A 
QYMO1A 
QYMO1B 
QYMO1B 
QYMO1B 
QZNO 1B 
QZNO 1B 
QZNO1B 
QZNO1C 
QZNO1C 
QZNO1C 
QYNO1B 
QYNO1B 
QYNO1B 
QYNO1C 
QYNO1C 
QYNO1C 
QZO0 1A 


22 837E-08 
8.231E-08 
52071E-08 
3.298E-08 
7. 709E-08 
3. 664E-08 
2. 815E-08 
1.051E-07 
4.062E-08 
3.477E-08 
8.072E-08 
4.177E-08 
2 9I30E-08 
1. 301E-07 
52 590E-08 
3. 964E-08 
1. 583E-07 
4. 936E-08 
3.418E-08 
7. 064E-08 
4.177E-08 
3-693E-08 
2.487E-08 
2.2311E-08 
1.074E-08 
1. 191E-08 
Te 525E-09 
4e494E-09 
6. 7/06E-09 
4. 147E-08 
3. 298E-08 
1- 858E-08 
6. 230E-08 
3. 296E-08 
2 386E-08 
4. 624E-08 
1. O64E-08 
1. 607E-08 
52 286E-08 
2e 684E-08 
1. 884E-08 
52 IS4E-08 
2e 664E-08 
1.905E-08 
6. 815E-08 
3. 344E~-08 
2e219E-08 
4, 801E-08 
4.008E-08 
2e 724E-08 
6. 637E-08 
3. 381E-08 
2.367E-08 
9. 948E-08 


52.9 
51.8 
66.5 
57.8 
68.5 
72.0 
70.6 
Zits 
48.1 
51.7 
53.2 
53.8 
49.4 
3%. 2 
oe Pam 
53. 4 

924 
41.0 
48.1 
66.0 
40.3 
49.5 
56.3 
52.8 
34.3 
42.1 
37.6 
37.5 
63.4 
2629 
31.5 
32.0 
42.4 
53.0 
52.6 
9.2 
12.3 
12.8 
36.0 
39.8 
38.2.6 
29e8 
34.3 
fe as Pa 
41.4 
53.0 
49.7 
23.90 
25-2 
2504 
41.4 
48.0 
5325 
41.3 
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QZOO1A 
QZO01A 
QZ001B 
QZ00 1B 
QZ001B 
QYOO 1A 
QYOO1A 
QYOO 1A 
QYO0O1B 
QYO01B 
QYO0O1B 
QOZPO1A 
QZPO1A 
QZPO1A 
QZP01B 
QZP01B 
QZP01B 
QZPO01B 
QZPO1B 
QYPO1B 
QYPO1B 
QYPO1B 
~ QYPOIC 
QYPO1C 
QYPO1C 
QYPO1C 
QZQO1A 
QZQ0 1A 
QZQO1A 
QZQ01B 
QZ001B 
QZ 00 1B 
QYQO1TA 
QYQOTA 
QYQO1A 
QYQO 1B 
QYQO 1B 
QYQO 1B 
QYQO1B 
QZRO1A 
QZRO1A 
QZRO 1A 
QZRO1B 
QZRO 1B 
QZRO1B 
QZRO 2B 
QZRO2B 
QZRO2B 
QZRO2D 
QZRO2D 
QZRO2D 
QYRO1A 
QYROTA 
QYROMA 


13.43 
13.43 
13.77 
13.77 
13.77 
13.23 
13.23 
13.23 
13.95 
13.95 
13.95 
13. 48 
13.13 
13.43 
13.10 
13. 10 
13.10 
13.10 
13. 10 
13.40 
13.40 
13.40 
13.54 
13.54 
13.54 
13.54 
13.34 
13. 34 
13.34 
13.53 
¥3~:53 
13.93 
13. 33 
135,33 
13. 33 
13.59 
13.59 
T3eoo 
13.59 
13. 38 
13. 38 
13.38 
13.39 
13.39 
13.39 
13. 31 
13.31 
13. 31 
13245 
13.45 
13.45 
13.90 
13.90 
13.90 


-150.0 
=~15020 
-150. 0 
- 150. 0 
-150.0 
—165.:0 
-165.0 
- 165.0 
-165.0 
-165.0 
-165. 0 
-180. 0 
-180.0 
-180. 0 
-180. 0 
-180.0 
-180. 0 
-180. 0 
-180.0 
199. 0 
Ged oe 8) 
= 9 54.0 
= 19 56i0 
oe er, 
-195. 0 
= 795.0 
-210. 0 
-210. 0 
-210.0 
-210.90 
-210.0 
21 GeO 
eye are) 
meee, 
22000 
2255.0 
Oe Dak 
—225..0 
=22 500 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
-240.0 
~259~.0 
-255. 9 
=O 5 se0 


4.914E-08 
32 1895-08 
1. 142E-07 
5. 534E-08 


. 3 739E-08 


1. 258E-07 
3. 230E-08 
1. 736E-08 
1.181E-07 
5.503E-08 
3. 923E-08 
5. 704E-08 
22 149-08 
1.916E-08 
12467E-07 
2.904E-08 
22014E-08 
1.782E-08 
1.093E-08 
8. 265E-08 
3. 3458-08 
2. 230E-08 
1.007E-07 
4.237E-08 
2. 582E-08 
1. 794E-08 
8. 445F-08 
2.678E-08 
1. 660E-08 
9.026E-08 
3. 584E-08 
2.422E-08 
8.819E-08 
2-737E-08 
2.543EF-08 
9.321E-08 
4. 156E-08 
3.201E-08 
2. 254E-08 
8. 861E-08 
3.537E-08 
22 258E-08 
9.497E-08 
2. 864E-08 
1.517E-08 
1. 009E-07 
4. 349EF-08 
2.457E-08 
9.673E-08 
3. 886E-08 
1.690EF-08 
8.793E-08 
1.953E-08 
1. 764E-08 


49.3 
5203 
3502 
53.4 
42./ 
301.4 
Oo2 
36.6 
18.6 
36.6 
37.0 
21.8 
54.0 
52.7 
20.2 
28.5 
586i2 
41.6 
63.9 
44,3 
60.5 
60.0 
35.9 
30.5 
55.0 
61.90 
324.3 
10.8 


. 36.28 


33.6 
56.6 
58.9 
327.9 
41.0 
58.4 
64.1 
69.8 
70.9 
sas 
4.3 
28.3 
36.1 
29.6 
47.1 
47.8 
3 oe a 
65. 3 
66.22 
34.6 
19.7 
7.2 
320.2 
24.7 
32a: 
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QYROI1A 13.90 -255.0 30.0 1.483E-08 38.6 23.1 
QYRO1B 13.40 -255.0 0.0 9.529E-08 31.4 17.9 
QYRO1B 13.40 -255.0 10.0 3.356E-08 39 2 16.3 
QYRO1B 13.40 -255.0 20.0 1. 460E-08 35.5 14.5 
QYRO1C 13.64 -255.0 0.20 7.2015E-08 18.6 40.7 
QOYRO1C 13.64 -255.0 10.0 22 ID0E-08 52.9 28.4 
QYRO1C 13.64 -255.0 20.0 1-2 751E-08 58.5 oe Pes 
QYRO1C 13.64 -255.0 30.0 9. 099E-09 56.8 36.0 
QYRO1C 13.64 -255.0 40.0 1.094E-08 359, 3 56.7 
QYRO1C 13.64 -255.0 60.0 6. 548E-09 47.2 19.2 
QYRO1C 13.64 -255.0 80.0 5e 198E-09 318.6 56.6 
QYRO1C 13.64 -255.0 100.0 8. 605E-09 72.3 31a 3 
QyYRO1C 13.64 -255.0 140.0 4.208E-10 314. 3 7022 
QYRO1C 13.64 -255.0 180.0 4.579E-09 8.8 5e.2 
QZS01B 13.40 -270.0 0.0 9. 943 E-08 36.5 15.4 
QZSO1B 13.40 -270.0 10.0 3. 862E-08 46.5 8.3 
QZSO1B 13.40 -270.0 20.0 2-036E-08 48.1 3.7 
QZS01C 13.30 -270.0 0.0 1.061E-07 bs 36.2 
QZSO1C 13.30 -270.0 10.0 4. 605E-08 69.8 40.7 
QzZSO01C 13.30 -270.0 20.0 3.842E-08 60.9 89.13 


Choteau Section 


QUFO1B 12.86 11.0 0.0 6. 523E-08 70.4 63.0 
QUFO1B 12.86 11.0 5.2 0 5. 865E-08 80.4 67.2 
QUFO1B 12.86 11.0 10.0 5. 237E-08 81.4 67.7 
QUFO1B 12.86 11.0 20.0 4. 536E-08 104.8 25g 
QUFOIC 13.13 11.0 0.0 6. 242E-08 78.9 66.2 
QUFOI1C 13.13 11.0 Zuid 52 17I9E-08 89.4 68.8 
QUFOIG . ¥3.33 11. 0 520 5. 436E-08 89.7 7122 
QUFO1C 13.13 11.0 10.0 4.994 E-08 95.5 69.6 
QUFO1C 13.13 11. 0 20.0 4.162E-08 97.6 68.3 
QUFO1C 13.13 11.0 30.0 32471E-08 108.6 69.6 
QTUOIA 12.29 9.0 0.0 6. 752E-08 12.8 51.0 
QTUOTA 12.29 9.0 265 6. 802E-08 9.7 50.7 
QTUOIA 12.29 9.0 520 6.570E-08 8.9 49.8 
QTUOIA 12.29 9.0 10.0 6. 081E-08 8.0 49.1 
QTUOIA 12.29 9.0 20.0 5. 169E-08 7.0 47.9 
QTUOIA 12.29 9.0 30.0 4, O34E-08 Se 1 47.7 
QUEOIA 12.71 +00 0.0 5. 466E-08 Zee 59.4 
QUEOIA 12.71 -7.0 529 5.044E-08 24.0 62.9 
QUEOQ1A 12.71 -7.0 10.0 4.617E-08 21.8 65.4 
QUEO1A 12.71 +TeQ 20.0 3.821E-08 27.9 67.3 
QUEO1B 12.64 -7.0 0.0 5-491E-08 23.6 60.8 
QUEO1B 12.64 She 0 2a5 52 134E-08 Z2el 64.7 
QUEO1B 12.64 -7.0 5.0 4.896E-08 20. 6 65.1 
QUEO1B 12.64 -7.0 10.0 4.498E-08 17.3 65. 6 
QUEO1B 12.64 -7.0 20.0 3. 6545-08 24. 2 66.5 
QUEO1B 12.64 =750 30.0 3. 087E-08 1957 66. 6 
QUDO1IA 11.92 -17.0 0.0 4.400E-08 -15.0 75.6 
QUDO1A 11.92 -17.0 10.0 4. 224E-08 -37.4 2303 
QUDO1A 11.92 -17.0 20.0 4. 288E-08 -42.8 74.5 
OUDO1B 3.72 --72.0 0.0 6.723E-08 19.4 Bes Be 
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QUDO01B 13.272 -17.0 
QUDO01B 13.72 -17.0 
QTAOIA 11.83 -20.0 
QTAOTA 11.83 -20.0 
QTAOIA 11.83 -20.0 
OTAO1B 11.46 -20.0 
QTAOQ1B 11.46 -20.0 
QTAO1B 11.46 ~-20.0 


= 


5. 809E-08 = 85.1 
5. 886E-08 212.6 81.4 
6. 367E-11 343.0 50.0 
52 222E-08 344.7 50.4 
4.481E-08 344.2 49.0 
6.308E-08 342.2 49.7 
5e 423 E-08 340.4 50.3 
4. 445E-08 338.4 49,3 


N 


©Ooqgqoooooceqooorocoococnc”eqoqooqoncoc$coododooqoooqoooqonooocncomooqooqooqoqoqoonoeoq0qoceooano0o°cd ea © 


tN) = 
e 6 


od 


N 


QUCO1B 13.47 -30.5 ° 6. 380E-08 ~42.5 Fe 
QUCO1B 13.47 -30.5 1 4. 899E-08 =a 65.9 
guco1B 13.47 -30.5 2 3.687E-08 =€05/5 67.0 
OUGE@TE 13.80 ~=30.:5 6. 846E-08 ~35.9 T1207 
QUCOIE 13.50 *=30.5 1 5.601E-08 -48.7 63.6 
Qgucoic 13.50 -30.5 2 3.910E-08 -46.3 60. 3 


52 1215-08 347.9 57.7 
5. 304E-08 344.4 5 8. 3 
4. 961E-08 342.8 57.9 
4.707E-08 327.6 47.6 
5= 532E-08 0.3 60.2 
5. 433E-08 357.7 60.8 
5. 360E-08 356.4 61.7 


QTTO1B 12.09 -42.0 
QTTO1B 12.09 -42.0 
QTTO1B 12.09 -42.0 
QTTO1B 12.09 -42.0 
QTTO1C 12.48 -42.0 
QTTO1C 12.48 -42.0 
QTTO1C 12.48 -42.0 


ND = 


QTTO1C 12.48 -42.0 1 4. 812E-08 357.0 61.6 
QTTOIC 12.48 -42.0 2 4.076E-08 354.9 59.0 
OTTOIC - 12.498 #=-42.0 3 3.471E-08 357.4 Pe 
QUBOITA 13.07 -47.0 5e 101E-08 =5 363 64.3 
QUBO1A 13.07 ~-47.0 1 4. 330E-08 ~52Z.9 62.6 
QUBO1A 13.07 -47.0 z 4.286E-08 -54.6 58.8 
QUBO1B 12.20 -47.0 5. 8S45E-08 -40.6 D5a4 
QUBO1B 12.20 -47.0 1 4. 751E-08 =52. 7 66.9 
QUBO1B 12.20 -47.0 2 3.412E-08 -56.9 59.4 
QUAO1A 13.51 -6€0.0 5e 563E-08 -0.5 74.0 
QUAOTA 13.51 -60.0 1 5. 110E-08 =Zoao 7422 
QUAOQ1A 13.51 -60.0 Zz 4,.340E-08 ~S1eg 7729 
QUA01B 13.81 -60.0 4.012E-08 -4,1 5963 
QUA01B 13.81 -60.0 1 4, 762E-08 =250§ 58.1 
QUAO01B 13.81 -60.0 2 52 751E-08 -~34.0 59.7 


5.2 618E-08 347.3 62.3 
4. 769E-08 341.8 64.8 
4.029E-08 333.0 65.1 
5< 227E-08 348.3 52e5 
4.400E-08 SS ta) 54.8 
3.598E-08 336.7 54.3 
6.630E-08 337.1 50.0 
6.053E-08 332.4 52.6 
5.-615E-08 326.4 52.1 
4,678E-08 327. 2 Daa 


QTSO1A 12.43 -62.0 
QTSO1A 12.43 -62.0 
QTSO1IA 12.43 -62.0 
OTSOIB, 13.33 *=62.0 
OFSOTB 13.33 ~=-62..0 
QOTSOIB 13.33 ~=62.0 
QTQOIA 12.43 -78.0 
QTOOIA 12.43 -78.0 
QTQOIA 12.43 -78.0 
QTQO1A 12.43 -78.0 


ome ht — 


th) = 


QOTQOIB 12.53 -=78.0 6. 364E-08 —Viad 43.7 
QTQQ01B 12.53 -78.0 1 5. 862E-08 -26.6 48. 8 


QTROIA 12.53 -106.0 
QTROIA 12.53 -106.0 
QTROIA 12.53 -106.0 
QTROIA 12.53 -106.0 
QTRO1B 12.05 -106.0 


5. 486E-08 3595.4 51. 1 
4. 883E-08 356.8 522 
4.492E-08 35849 52.4 
3.895E-08 358.0 49.28 
4.486E-08 33362 554.9 
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QTRO1B 
QTRO1B 
QTRO1B 
OTRO1B 
QTRO1B 
OTOB6A 
QTOB6A 
OTOB6A 
QTOB6A 
QTOB6A 
QTOB6B 
QTOB6B 
QTOB6B 
QTOB6B 
OTOB6B 
QTOBSA 
QTOB5A 
QTOB5A 
OTOB5B 
QTOB5B 
OTOBSB 
QTOBYA 
OTOB4A 
OTOB4SA 
QTOB4B 
QTOBYB 
QTOBYB 
QTOB2A 
QTOB2A 
OQTOB2A 
QTOB2A 
OTOB2A 
QTOB2B 
QTOB2B 
QTOB2B 
QTOB2B 
QTOB2B 
QTOB2C 
QOTOB2C 
QTOB2C 
QTOB2C 
QTOB2C 
QTOB1B 
QTOB1B 
QTOB1B 
QTOB1C 
QTOBIC 
OTOB1C 
QTPO1B 
QTPO 1B 
QTPO1B 
OTP01B 
QTPO1C 
QTPO1C 


12.05 
12.05 
12.05 
12.05 
12.05 
ua 95 
13.895 
T3695 
13.95 
13.95 
12.42 
12.42 
12.42 
12.42 
12.42 
11.83 
11.83 
11. 83 
17.54 
11.54 
11.54 
12.08 
12.08 
12.08 
12.15 
12.15 
12.15 
12.94 
12.94 
12.94 
12.94 
12.94 
13.07 
13.07 
13.07 
13.07 
13.07 
13. 76 
13.76 
13.76 
13.76 
13.76 
10.43 
10.43 
10.43 
10.71 
10.71 
10.71 
10.79 
10.79 
10.79 
10.79 
11.61 
11.61 


-106.0 
-106. 0 
-106.0 
-106.0 
-106.0 
-120.0 
—T29.2,0 
-129.0 
-129.0 
-129.0 
~129.9 
-129.:9 
-129.90 
-129.90 
-129. 0 
-145.0 
-145.0 
-145.0 
-145.0 
-145.0 
-145.0 
-7160. 0 
- 160.0 
-160.0 


-160.0 


-206. 0 
-206.0 
-206.0 
-206.0 
-206.0 
-206.0 
-206.0 
-206.0 
-206. 0 
-206.0 
-206. 0 
-206.0 
-206.0 
-206.0 
-218. 0 
-218.0 
-218.0 
-218.0 
-218.0 
-218.0 
-250.0 
-250.0 
-250. 0 
-250.0 
-250.0 
-250.0 
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68.5 
67.8 
70.9 
76.1 
69.9 
66.7 
72.9 
80.4 
60.4 
63.1 
62.6 
60.7 
69.6 
68.3 
66.6 
47.0 
55.1 
5546 
54.7 
61.2 
63.0 
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APPENDIX 3 
The data used in Chapter 3 is listed here in the same 
format used Appendix 2. Distances are measured from the 


Cherryville Tephra. 


Riggins Road Data 


QJA01 ie 0.0 0.0 1,.401E-06 8.1 54.3 
QJA01 9692 0.0 205 12335E-06 9. 1 5309 
QJA01 Gene 0.0 Del 1.2237E-06 11.2 Saad 
QJA01 9252 0.0 10.0 1.019E-06 12.2 47.9 
QJA01 9252 0.0 15.0 9.230E-07 15.0 49.3 
QJA01 0092 0.0 20.0 7.622E-07 16.6 47.0 
QJA01 9.252 0.0 30.0 1.520E-06 18.2 T2eaf 
QJAO1 9532 0.0 40.0 203 16E-O7 22.7 44.4 
QJA02 9.68 0. 0 0.0 1.624E-06 1.7 63.2 
- QOJA0D2 9.68 0.0 10.0 1.245E-06 4.6 59.20 
QJBO1 10.03 -15.0 0.0 3. 649E-07 7.1 64.0 
QJBO1 10.03 =75.0 10.0 3- 181E-07 11.4 61.4 
QJBO2 10.91 15.0 0.0 3.2 940E-07 14.2 66.8 
QJBO2 10.91 z15e0 10.0 3.2363E-07 152.2 64.1 
QJCO1 11.44 -30.0 0.0 1.4 74E-07 356.5 65.3 
QJCO1 11.44 -30.0 10.0 1.248E-07 0.2 63.3 
QJCO2 10.93 -30.0 0.90 2- 103E-07 aek 699 
QJCO2 10.93 -30.0 10.0 1.484E-07 9.0 66.7 
QJDO1 10.68 -46.0 0.0 22 521E-06 8.5 66.9 
QJD0O01 10.68 -46.0 10.0 2-356E-06 8.0 66.6 
QJDO02 10.78 -46.0 0.0 2.2 261E-06 7.3 69.0 
QJDO2 10.78 -46.0 10.0 2. 133E-06 7.8 68.8 
QJEO1T 10.43 -64.0 0.0 9. 361E-07 721 7324 
QJEO1 10.43 -64.90 10.0 Se 059-07 6.8 7361 
QJEO2 10.34 -64.0 0.0 7.2209E-07 14.2 68.6 
QJEO2 10.34 -64.0 10.0 6. 989E-07 13.8 56.2 
QJFO1T 10,44 -88.0 0.0 2.047E-07 3.4 63.0 
QJFO1 10.44 -88.0 10.0 1.339E-07 2-4 5924 
QJFO2 10.46 -88.90 0.0 1.895E-07 1.5 71.1 
QJFO2 10.46 -88.0 10.0 1. 151E-07 0.9 68. 4 
QJGO1 10.23 -108.0 0.0 4. 092E-07 Sane 71.4 
QJGO1 10.23 -108.0 10.0 3.268E-07 4.8 712.2 
QJGO2 10.60 -108.0 0.0 32 244E~07 2.1 67.4 
QJGO2 10.60 -108.0 10.0 22 987E-07 6.3 67.0 
QJHO1 9232 —122.0 0.0 7. 362E-07 Se2 70.8 
QJHO1 9632 ~-722..0 2.065 Te 156E-07 726 70.4 
QJHO1 9.322 -122.0 5.0 6. 940E-07 725 70.4 
QJHO1 Gaas  ~b22..0 10.0 6.4 77E-07 7.9 69.7 
QJHO1 S632 —-12240 15.0 5 162E-07 7.9 69.4 
QJHO1 Gaga “2250 20.0 5- 183E-07 7.8 69.1 
QJHO1 9.32 —-322..0 30.0 3.739E-07 9.3 66.8 
QJHO1 9252 ~122..0 40.0 22329E-07 10.1 70.9 
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QJHO2 10.77 -—122.,0 0.0 1.4 74E-06 12.8 68.8 
QJHO2 10.77 -122.0 10.0 1. 184E-06 To? 64.9 
QJI01 10219 -137.0 0.0 8. 113E-07 Ss 68.1 
QJI01 10.19 -137.0 10.0 6.618E-07 0.9 66.9 
QJT0O2°.10.35 -337.0 0.0 3.308E-07 6.6 66. 6 
QJIO2 10.35 -137.0 10.0 2e 681E-07 4.0 64.1 
QJJ01 10.15 -155.0 0.0 Ze 255E-06 353.0 74.8 
QJJ01 10.15 -155.0 10.0 1. 825E-06 3525'6 7204 
QJJO2 10.21 -155.0 0.0 1. 407E-06 357.6 70.9 
QJJO02 10.21 -155.0 10.0 1. T15E-06 360.0 68.2 
QJKO1 10.00 -171.0 0.0 Je 333E-U6 7.26 7105 
QJKO1 10.00 -171.0 10.0 1.980E-06 tad 70.8 
QJKO2 9.92 - 11.0 0.0 1. 894E-06 1332 70.1 
QIKO2 9292 ~—171.0 10.0 1. 584E-06 14.8 69.8 
QJLO1 9.89 —“—185.90 0.0 6. 314E-07 12.3 69.2 
QJLO1 5.68 ~26520 10.0 4.286E-07 13.6 67.8 
QJLO2 10.00 -188.0 0.0 6. 365E-07 6.1 67.0 
QJLO2 10.00 -188.0 10.0 4.305E-07 7.5 64.0 
QJMO1 10.26 -208.0 0.0 1. 341E-06 3.6 72.3 
QJM01 10.26 -208.0 10.0 1. 045E-06 0.5 70.7 
QJMO02 10.40 -208.0 0.0 9..938E-C07 347.5 76.7 
-QGIMO02 10.40 -208.0 10.0 72937E-07 35123 W502 
QJNO1 10.39 -224.0 0.0 1.449E-06 1.1 69.8 
QJNO1 10.39 -224.0 10.0 1.233E-06 359.6 68.6 
QJNO2 10.32 -224.0 0.0 2- 263E-06 SD a2 71.6 
QJNO2 10.32 -224.0 10.0 1. 956E-06 92 69.1 
QJ001 7285 -—244.0 0.0 1.652E-06 4.1 7162 
QJ001 7.85 -244.0 10.0 1. 389E-06 0.2 69.6 
QJO02 9255 -244.0 0.0 1.817E-06 12.9 71.8 
QJO02 9.55 -244.0 10.0 1. 482E-06 Ded 68.4 
QJPO01 10.96 -268.0 0.0 2e499E-07 1.1 70.9 
QJPO01 10.96 -268.0 2.5 2e330E-07 2-6 7124 
QJPO01 10.96 -268.0 5.90 2e 201E-07 1.7 70.7 
QJPO1 10.96 -268.0 10.0 20 137E-07 3592.3 6359 
QJPO1 10.96 -268.0 15.0 1. 924E-07 0.6 69.4 
QJPO01 10.96 -268.0 20.0 1.631E-07 357.2 70.24 
QJPO1 10.96 -268.0 30.0 1.172E-07 Sa 66.2 
QJP01 10.96 -268.0 40.0 1.087E-07 28.9 41.1 
QJPO2 11.02 -268.0 0.0 4,.128E-07 0.1 70.6 
QJPO2 11.02 -268.0 10.0 3.437E-07 357. 1 70.3 
QJQ01 10.36 -289.0 0.0 3.255E-07 344.9 66.9 
QJQ01 10.36 -289.0 10.0 22465E-07 347.0 64.7 
QJQ02 10.48 -289.0 0.0 1.682E-07 3e5 70.0 
QJQ02 10.48 -289.0 10.0 1.180E-07 TeT 67.9 
QJRO1 10.91 -304.0 0.0 2.409E-07 7.8 72.3 
QJRO1 10.91 -304.0 10.0 1. 423E-07 6.9 71.1 
QJRO2 11.00 -304.0 0.0 20 244E-07 14.3 74.2 
QJRO2 11.00 -304.0 10.0 124 72E-07 13.2 72.1 
QJSO1 10.245 -328.0 0.0 1.425E-907 20.9 75.0 
QJSO1 10.45 -328.0 10.0 9.293E-08 14.9 70.6 
QJSO2 10.41 -328.0 0.0 1.671E-07 363 B29 
QJSO2 10.41 -328.0 10.0 8. 585E-08 Ze9 POLS 
QJTO1 10.67 -346.0 0.0 22331E-07 11.9 67.5 
QJTO1 10.67 -346.0 225 2e 125E-07 11.6 66. 8 
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QJTO1 10.67 -346.0 529 22 094E-07 11. 1 65.4 
QJTO1 10.67 -346.0 10.0 1.618E-07 929 66.9 
QJTO1 10.67 -346.0 15.0 1.036E-07 9.6 66.2 
QJTO1 10.67 -346.0 20.0 6.219E-08 11. 4 64.0 
QJTO01 10.67 -346.0 30.0 4 .946E-08 10.4 65.3 
QJTO1 10.67 -346.0 40.0 30 127E-08 9.0 69.7 
QJTO2 10.56 -346.0 0.0 Ze Q33E-O07 6.8 68.4 
QJTO2 10.56 -346.0 10.0 1.362E-07 ane 65.5 
QJUO1 10.40 -365.0 0.0 2- 312E-07 358.0 75.1 
QJU01 10.40 -365.0 10.0 1.622E-07 357. 6 74,1 
QJU02 10.76 -365.0 0.0 1.972E-07 0.1 64.4 
QJU02 10.76 -365.0 10.0 1.280E-07 358.9 61.8 
QIKO1 10.44 -380.0 0.0 22 9013E-07 14.1 66.5 
QIKO1 10.44 -380.0 10.0 1.610E-07 11.2 64.1 
QIKO2 10.66 -380.0 Q.0 1.812E-07 3.1 1249 
QIKO2 10.66 -380.0 10.0 1.315E-07 3.0 67.6 
QJVO1 10.88 -400.0 0.0 1.405E-07 359; 6 70.0 
QJV01 10.88 -400.0 10.0 8.296E-08 359.0 69.5 
QJVO2 10.42 -400.0 0.0 1.231E-07 357.4 59.6 
QJVO2 10.42 -400.0 10.0 Te O41E-08 354.6 ate! 
QJWO01 9.73 ~#19..0 0.0 3.167E-07 357.4 66.9 
~ QIW01 9.73 —419.0 2.5 2 914E-07 357.1 65.0 
QJW01 9.73 -419.0 5.20 2~e 639E-07 32369 63.4 
Quwo1 $2.73 -—499..9 10.0 2-011E-07 354.4 62.1 


QJWO1 59.73 -~4&19.0 15.0 1.415E-07 356.7 5h. F 
QJUW01 9.73 : ~419.0 30.0 6.2 O084E-08 336.5 69.5 


QJuw01 9.73 -419.0 40.0 4..063E-08 1. 1 48.1 
QJWO2 e971 —4499.0 0.0 2e380E-07 13.0 67.8 
QIWO2 9.97 -419.0 10.0 1.210E-07 6.3 61.7 
QJX01 9.56 —-439.0 0.0 6.518E-07 3.4 63.1 
QJX01 9,96 —§39.0 10.0 4.113E-07 204 60.9 
QJX02 10.03 -439.0 0.0 6.427E-07 8.5 65.1 
QJXO02 10.03 -439.0 10.0 4.031E-07 rs) 63.1 
QJY0O1 9.260 -459.0 0.0 1.347E-07 17.0 63.9 
QJYO01 9.60 -459.0 10.0 8.951E-08 15.5 57.8 
QJYO2 8.10 -459.0 0.0 1.226E-07 15.6 77.9 
QJYO2 8.10 -459.0 10.0 6.637E-08 13.9 74.3 
QJZ01 3.99 -502. 0 0.0 72 572E-08 18.9 63.1 
QJZ01 8.99 -502. 0 2.25 6. 983E-08 16.4 61.2 
QJZ01 6.99) —5025:0 5.0 6.209E-08 12.0 61.2 
QJZ01 9299  ~=-502./0 10.0 4. 854E-08 17.3 59.0 
QJZ01 8.97. -502.70 15.0 4.214E-08 13.6 54.0 
QJZ01 8699 —~502.0 20.0 3.461E-08 22.0 75.0 
QJzZ01 8.99, —502. 0 30.0 22 S24 E-08 19.1 56.6 
QJZ01 §.99 -502.0 40.0 1.817E-08 71.0 79.1 
QJZ02 8.77 -502.0 0.0 4.532E-08 7.7 65.9 


QJZ02 S237 =502..0 10.0 22 552E-08 35 t..2 60.7 
QIA0O1 10.69 -519.0 0.0 1. 148E-07 16.0 59.6 
OQIAOT 10.69 ~5179.0 10.0 6./782E-08 13.2 57.0 
QIAO02 10.98 -519.0 0.0 9. 545E-08 21.1 5 8.0 
QIAO2 10.98 -519.0 10.0 4. 273E-08 A2ad 57.3 
QIBO1 10.24 -540. 0 0.0 22899E-07 12.6 7420 
0.0 


QIBO1 10.24 -540.0 10, 1.401E-07 11.9 644.3 
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QIBO2 10.33 -540. 0 2. 215E-07 Zane 78.3 
QIBO2 10.33 -540.0 1 9. 641E-08 15. 6 71.7 
QIco1 10.04 -557.0 8.892E-07 6.7 66. 4 
QICO1 10.04 -557.0 1 6.968E-07 6.7 63.3 
QICO2 9.94 -557. 0 1. 536E-06 6.8 71.6 
QIC02 969% 55720 1 1.274E-06 4.4 70.1 
QIDO01 10.87 ~576.0 52 210E-07 559 71.5 
QIDO01 10.87 -576.0 4. 346E-07 1. 6 70.0 


OID@2 10545 -S76.0 
QIDO02 10.45 -576.0 
QIEO1 11.38 -596.0 
QTEO1 11.38 -596.0 


22288E-07 358.5 69.8 
1. 820E-07 349.1 67.9 
9.711E-08 358.6 TGs9 
7. 281E-08 356.3 1523 
QIEO1 11.38 -596. 0 6. 975E-08 35745 74.2 
QIEO1 11.38 -596.0 1 59 99E-08 3543.2 713.7 
QTEQJ 41.238 —596.0 15.0 4.902E-08 a5 580 79.4 
QIEO1 11.38 -596.0 20.0 5. 9035-08 43.6 56. 8 
OLHOT 11.38 + —59650 30.0 52 967 E-08 27.7 17.7 
QIEO1 11.38 -596.0 40.0 22629E-08 353.0 60.1 
QIEO1 11.38 -596.0 50.0 1.659E-08 25.6 69.6 

QIEOT 11.38 -—596. 0 60.0 4.380E-09 306.2 ios, 


—_ = 
OmMmNCeCOCOCOCOCoOOCoO eo 
oCoomooqnoeoeococneeoo 


QIEO2 10.71 -596.0 0.0 22 249E-07 6.6 7724 
‘OQLE02 10071 )5-596..0 10.0 1.503E-07 0.6 74.6 
OTEO2> 10577 ~=5962:0 15.0 1. 240E-07 1.8 73.0 
QIEO2 10.71 -596.0 17.65 1.025E~07 4.0 71.7 
QIEO2 10.71 -596.0 20.0 Veto e-O7 929 aon) 
OTEOZ2 10.71 --596.0 2205 1-027E-07 26.25 72.7 
OPEOZe? 1025/10 = 5965.0 25.0 953032 -08 6.5 73.26 
QOLEO2 (10571 —-59650 Ziad 9. 083E-08 15.5 7329 
QIEO2 10.71 -596.0 30.0 8.779E-08 10.3 71.4 
QIEO2 10.71 -596.0 32.5 8. 195E-08 11.5 72.7 
QIEO2 (105717 -596. 0 35.0 72 190E-08 357. 3 71.6 
QIEO2 10.71 -596.0 3725 7. 100E-08 63.1 65.2 
OLEQ2 10571 ~-396..0 40.0 6.521E-08 Zoo 69.7 
QIEO2 10.71 -596.0 42.5 6. 949E-08 1. 6 Bath 
QIEO2 10.71 -596.0 45.0 5. 760E-08 1.8 69.4 
QIEO2 10.71 -596.0 47.25 5e 292E-08 13.7 To0 | 
QIEO2 10,71 -596..0 50.0 4.550E-08 dod 70. 3 
QOLEO2 10571 -596.0 52.25 4. 835E-08 6. 6 71.4 
OLEO2Z 9 A021) — 5965.0 55.0 32507E-08 14.3 70.2 
OLTEO2. 102717) —3596. 0 60.0 20093808 Zee 67. 1 
QIEO2 10.71 -596.0 65.0 2e270E-08 341.6 66.5 
QIEO2 10.71 -596.0 70.0 1.384E-08 39.3 T3.9 
QIFO1 10.98 -616.0 0.0 1.181E-07 10.0 75.4 
QIFO1 10.98 -616.0 10.0 7a 326E-08 6. 4 71.8 
OLFO2) T1356, —6:16..0 0.0 1.068E-07 1.0 80.9 
OLFO2 T1256 “-616. 0 10.0 6. 841E-08 4.5 73.6 
QIGO1 10.51 -634.0 0.0 1.2776E-07 15.0 3a) 
Q1GO1 10.51 -634.0 10.0 1. 134E-07 12.7 70.0 
QIGO2 10.64 -634.0 0.0 12 7195-0T 20.9 72.0 
QIGO2 10.64 -634.0 10.0 12115E-07 14.1 68.6 
QOTHO?T 10586 -655.0 0.0 2a 622E—07 3.2 75.0 
QIHO1 10.86 -655.0 225 22 240E-07 Se 1 74,9 
QIHO1 10.86 -655.0 10.0 1.681E-07 6.6 73.) 
QITHO1 10.86 -655.0 15.0 1.347E-07 8.0 68.7 
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QIHO1 10.86 -655.0 20.0 1.058E-07 Gere 126 
QIHO1 10.86 -655.0 30.0 72916E-08 14,7 48.1 
QTHO?, 10.868°=65550 40.0 5. 204E-08 10.2 5S eio 
OQTHO2 10.3398+65520 0.0 1. 149E-07 333.2 83.29 
QOPHO2- 10533) =—6552.0 10.0 De Z91E-08 8.9 7725 
QLHO2ZSPI02338% (26555.0 15.0 4.047E- 08 358. 8 7822 
QITHO2 10.33 -655.0 20.0 2. 672E-08 6.1 76.5 
QII01 9.83 -670.0 0.0 1.348E-07 7.8 70.3 
QOiI01 9.83 -670.0 10.0 7. 203E-08 Ze l 65.4 
QITO2 410081) -67020 0.90 22005E-07 Ts2 74.0 
QIIO02 10.81 -670.0 10.0 1. 0835-07 3.8 iZare 
QIJ071 11.01 -686.0 0.0 Te O46E-07 3e9 7146 
QTJ01 11.01 -686. 0 10.0 4,962E-07 1.8 68.9 
QIJ02 11.28 -686.90 0.0 82460E-07 8.3 62.7 
QIJ02 °11.28 -686.0 10.0 52 971E-07 5ef 60.8 
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APPENDIX 4 


FORTRAN Listing of Synthetic Data 
Generating Program "WRTFREQ" 


COMPLEX CX (50) ,CX1,CX2,CX3 

This program (WRTFREQ) produces a synthetic sequence of 
paleomagnetic directions (specified length N) 
containing a maximum of three superimposed looping 
Signals of arbitrary frequency, sense, 

ellipticity, and amplitude. 

Fisherian noise of arbitrary level may 

be added to the data. 

The data can be produced at non uniform intervals of 
time if a file containing N times is provided through 
unit 2 (format (7£9.4)). 

Output is onto a file attached to unit 1 andisa 
sequence of declinations and inclinations suitable for 
input into Maxent (2£9.4). 

If this program is compiled along with the attached 
subroutine (Rotate) then it can be run as follows: 


R -LOAD#+*IMSLLIB 1=OUTPUT 2=TIMES 


REAL LEVEL (50) ,INC,DINC2 (50) , THETA (50) 
SCALE=1.0 

RAD=57. 29576 
PT=2.0*ARCOS (-1. 0) 
WRITE (6,101) 

READ(5,100) P1,P2,P3 
WRITE (3, 102) 
READ(8,100) B1,B2,B3,A1,A2,A3,SD 
SD=SD/RAD 

SD= (1. -COS (SD) ) 
WRITE (6, 103) 
READ(5,100) AVEDEC,AVEINC 
AVEDEC=AVEDEC/RAD 
AVEINC=(90.-AVEINC) /RAD 
F1=1./P1 

F2=1./P2 

F3=1./P3 

WRITE (6, 105) 

READ(7,106) IX 

CALL GGUB(IX,1,VV) 
W1=F1*PI 

W2=F2*PI 

W3=F3*PI 

WRITE (6, 104) 

READ(5,100) SPACE 

IF (SPACE.EQ.0.0) GOTO 3 
DO 4 I=1,50 

READ(2,100) LEVEL(TI) 
SCALE=50.0/LEVEL (50) 
CONTINUE 

XTOT=0. 

CALL GGUB(IX,50, THETA) 
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CALL GGEXP (IX,SD,50,DINC2) 
DO 1 I=1,50 
IF(SPACE.EQ.0.0) LEVEL (I) =FLOAT (I) 
LEVEL (I) =LEVEL (1) *SCALE 
CX 1=CMPLX (COS (W1* LEVEL (I)) ,SIN(W1*LEVEL(I))) 
CX 2=CMPLX (COS (W2*LEVEL (I) ) , SIN (W2*LEVEL (I) )) 
CX3=CMPLX (COS (W3*LEVEL (I)) , SIN (W3* LEVEL (1) ) ) 
CX (I) =A1*CX14#B 1* CONJG (CX 1) #A2*CX2+ B2*CONJG (CX2) 
#4+A3*CX34+B3*CONJG (CX3) 
DEC=ATAN2 (-AIMAG (CX (I) ) ,REAL(CX(I))) 
INC= (90.-CABS (CX (I))) /RAD 
X=COS (DEC) *COS (INC) 
Y=SIN (DEC) *COS (INC) 
Z=SIN (INC) 
DINC1=PI/4.-INC 
DINC2 (I) =ARCOS (1.-ABS (DINC2(I))) +DINC1 
CALL ROTATE(X,Y,-DEC) 
CALL ROTATE(Z,X,-DINC2(T)) 
THETA (I) =THETA (I) *PI 
CALL ROTATE(X,Y, THETA (I) ) 
- CALL ROTATE (Z,X,DINC1) 
CALL ROTATE (X,Y, DEC) 
CALL ROTATE(Z,X,AVEINC) 
CALL ROTATE(X,Y,AVEDEC) 
DEC=ATAN2 (Y,X) *RAD 
INC=ARSIN (Z/SORT (X¥X+Y*Y+Z*Z) ) ¥*RAD 
1 WRITE(1,100) DEC,INC 
WRITE (4,106) IX 
WRITE(3,100) F1,F2,F3 
100 FORMAT(7F9.4) 
101 FORMAT(*ENTER 3 PERIODS (REAL) ' , 
102 FORMAT (*ENTER 3 CLOCKWISE THEN 3 COUNTERCLOCKWISE °, 
#*AMPLITUDES THEN STANDARD DEVIATION (DEGREES) ') 
103 FORMAT(" AVERAGE DECLINATION, AVERAGE INCLINATION;') 
104 FORMAT("ENTER 0.0 FOR EQUAL SPACING; °, 
#*ENTER 1.0 OTHERWISE') 
105 FORMAT(*ENTER IX‘) 


106 FORMAT (120) 


STOP 

END 

SUBROUTINE ROTATE(X,Y,ANGLE) 

SUBROUTINE ROTATE WILL ROTATE A UNIT VECTOR CLOCKWISE 
ABOUT THE X, Y, OR Z AXIS BY THE NUMBER OF RADIANS 
GIVEN IN THE INPUT PARAMETER "ANGLE". THE FIRST TWO 
PARAMETERS ARE THE VECTOR COMPONENTS WHICH WILL BE 
AFFECTED BY THE ROTATION. IT IS IMPORTANT THAT THEY 
BE ENTERED IN THE ORDER (Y FOLLOWS X, Z FOLLOWS Y, 
OR X FOLLOWS Z). THE MODIFIED COMPONENTS ARE STORED 
BACK INTO THESE TWO INPUT VARIABLES. 
XX=X*COS (ANGLE) -Y*SIN (ANGLE) 

Y=Y*COS (ANGLE) +X*SIN (ANGLE) 

X=XX 

RETURN 

END 
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Seed Integers Used in Chapter 4 


926686201 
748021105 
1448281358 
630620298 
840290746 
1680343038 
1310233148 
1397963139 
1154274150 
1116012542 
840290746 
1577741978 
1607589824 
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APPENDIX 5 


FORTRAN Listing of Complex Maximum Entropy Routine 
for Analysing Paleomagnetic Data 


COMPLEX CX (120),CB1(120) ,CB2(120) ,COM,CA (120), 
#CAA(120) ,CONJ 
This program (Maxent) calculates and plots the 
complex spectrum of a paleomagnetic time sequence. 
Input is through unit 1 and is assumed to be a file 
containing declination and inclination pairs 
(one pair to a line) in the format (2f9.2). 
Equal spacing in time between consecutive pairs 
is also assumed. 
The remanance directions are mapped onto 
the complex plane about their vector mean 
aS suggested by Denham (1975). 
The resulting time series is then analysed using the 
maximum entropy routine given in Anderson (1974) 
(Anderson's routine had to be modified slightly so 
that it would accept complex data). 
The program can Simply list the power densities (in 
degrees squared per unit frequency) associated 
with the principal peaks of the spectrum or it can 
produce a power density plot (on Calcomp) as well. 
The user may also integrate under spectrai peaks to 
obtain the power (amplitude squared) over a given 
frequency range. 
The FPE criterion of Ulrich and Bishop (1975) will 
be used to automatically choose a filter length 
appropriate for the data if the user wishes. 
If the automatic picking feature is used, then only 
filter lengths between 6 and N/2 (N = number of 
data points) will be considered. 
The FPE must also have a "proper" minimum within 
this range of filter lengths or no filter length 
will be chosen. 
If this program is compiled along with the attached 
subroutines it may then be run as follows: 
R ~LOAD#+* I MSLLIB+* PLOTLIB+ALAB: ARCHIVE 
+ALAB3O.BLOCKDATA 1=I 9=-P T=1M 
Where "I" is the input file of remanance directions 
and "-P" is the plot description file. 
REAL X (430) ,Y¥2 (430) ,Y (430) ,INC,Z (120) , F (50) ,P01(50) 
#, PO2 (3) ,PHS (50) 
INTEGER ITITLE (10) ,1R (430) 
DATA RAD,AVEX,AVEY,AVEZ,AVEDEC,AVEPOW,CC, XMIN, XMAX, 
#SCALE, ISPEC/57. 29576, 
#065 695-4 a7 05 700 fe 1EO7 F< 55007 1a, 07 
PT=6. 28 31853 
DO 1 I=1,100 
READ(1,99, END=10) DEC, INC 
DEC=DEC/RAD 
INC=INC/RAD 
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X (I) =COS {INC) *COS (DEC) 
Z (I) =SIN (INC) 

Y¥ (I) =COS (INC) *SIN (DEC) 
AV EX=AVEX+ X (I) 
AVEY=AVEY¢ Y (I) 
AVEZ=AVEZ+Z (I) 

N=I-1 

C=CC/FLOAT (N) 
AVEX=AINT (AVEX*C) /CC 
AVEY=AINT (AVEY*C) /CC 
AVEZ=AVEZ/FLOAT (N) 

IF (ABS (AVEX) +ABS (AVEY) .NE.0.0) AVEDEC= 


#-ATAN2 (AVEY, AVEX) 


AVEINC=-ARCOS (AVEZ/SQRT (AVEX**2+AVEY** 24+ AVEZ**Z2) ) 
DO 2 I=1,N 

CALL ROTATE (X(I) ,Y (1) , AVEDEC) 
CALL ROTATE (Z(1) ,X (12) ,AVEINC) 
DEC=ATAN2 (Y (1) ,X (1) ) 
INC=ARSIN (Z(1)) *RAD 

X (I) = (90.0-INC) *COS (DEC) 

Y (I) =(90.0-INC) *SIN (DEC) 

DO 3 I=1,N 

CX (1) =CMPLX (X (I) ,¥ (1)) 

DO 60 I=1,N 
AVEPOW=AVEPOW+CABS (CX (I) ) **2 
AVEPOW=AVEPOW/FLOAT (N) 
WRITE(6,107) AVEPOW 

WRITE (6,116) 

READ(5,101) IFPE 
IF(IFPE.EQ.0) GOTO 64 

MM=N/2 

GOTO 65 

WRITE (6, 102) 

READ(5,101) MM 

IF (MM.EQ.0) STOP 

nm=mam- 1 

WRITE (6, 108) 

READ(5,101) ISPEC 

CONTINUE 


COMMON CX,CAA,CB1,CB2,CA,COM,CONJ, Z,P0,PI,MM,N,ERR 


EXTERNAL POWER 

CALL MAXENT 

IF (ISPEC.EQ.1) GOTO 89 
IF(IFPE.EQ.0) GOTO 62 
MM=MM+1 

DO 61 I=1,MM 

IR(I)=1I 

IF(IFPE.EQ.1) GOTO 63 
IF(IFPE.GE.2) WRITE(6, 98) (IR(I) ,Z(1I) ,I=1,MM) 
IF(IFPE. EQ. 3) GOTO 63 
IFPE=0 

GOTO 64 

MN=MM-4 

CALL VSORTP (Z(5) ,MN,IR(5)) 
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KL=MM 

KM=5 

DO 49 I=5,MM 
IF(KM.GT.KL) GOTO 52 
IF(IR(I).EQ.KM) GOTO 51 
IF(IR(I).EQ.KL) GOTO 53 
GOTO 50. 

KM=KM+1 

GOTO 49 

KL=KL-1 

CONTINUE 

WRITE (6,150) 

FORMAT(*NO PROPER FPE MINIMUM CAN BE FOUND. PLEASE 


# PICK YOUR OWN FILTER LENGTH. ‘) 


STOP 

WRITE(6,119) IR(I) 
MM=IR(I)-1 

IFPE=0 

GOTO 66 

CB1(1)=(1.0,0.) 

MM2=MM+1 

DO 20 T=2,MM2 

CB1(1I)=-CA (I-1) 

CALL ZCPOLY (CB1,MM,CB2,1ER) 
DO 28 I=1,MM 

F (I) =ATAN2 (AIMAG (CB2(I)) , REAL (CB2 (I))) /PI 
PO1(1I) =POWER(F(I)) : 
COM=P0/COM 

PHS (I) =ATAN2 (AIMAG (COM) , REAL(COM) ) *RAD 
J=1I+360 

Y (J) =PO1(TI) 

X (J) =F (1) 

IR (J) =d 

DO 47 I=1,MM 

IR(1I) =I . 

CALL VSORTP(F,MM,IR) 
WRITE(6,120) (F (I) ,PO1(IR({I)) ,I=1,MM) 
WRITE (6,110) 

READ(5,101) IPLT 
IF(IPLT.EQ.0) STOP 

N=360+MM 

DO 18 'T=17360 

IR(I) =I 

X (I) =-2 5+. 5* FLOAT (I) /180.0 
Y (I) =POWER (X(1)) 

CALL USMNMX(Y,N,1,C,YMAX) 
CALL VSORTP(X,N,IR) 

DO 69 I=1,N 

Y2 (1)=Y¥ (IR (I)) 
TF(ERR.~EQ. 1.0) WRITE(6,103) 
DO'671I=1;776 
IF(YMAX.LT.10.**1) GOTO 68 
CONTINUE 

FOUR=4. *10** (I-2) 
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YMAX=FLOAT (INT (YMAX/FOUR) +1) * FOUR 
WRITE (6,114) 
READ(5,115) ITITLE 
CALL PLOTS 
CALL FACTOR (SCALE) 
CALL PLOT (1.0,1.0,-3) 
X (N+ 1) =XMIN 
X (N+2) = (XMAX-XMIN) /5. 
Y2 (N+1) =0.0 
Y2 (N42) =YMAX/4.0 
CALL FAXTS2 (0.0/0.0, * FREQUENCY *,-9,5.90,0.0,XMIN 
#,X (N-2) ,0.5) 
CALL AXIS2 (0.0,0.0,*POWER DENSITY',13,4.0,90.0 
#,0.0,Y2(N+2) ,-—1. 0) 
CALL LINE(X,Y2,N,1,0,003) 
CALL SYMBOL(1.0,4.1,.20,ITITLE,0. ,40) 
CALL PLOT(0.0,0.0, 999) 
STOP): 
89 WRITE (6, 106) 
4 READ(5,113) F1,F2,FINC 
TP (F1.EQ.F2) STOP 
NX=INT ((F2-F1) /FINC) +1 
FFI=F1 
TOT=0.0 
DO 31 J=1,NX 
TOT=TOT+GLO96 (POWER ,FF1,FF1+FINC,144) 
31 FFI=FFItFINC 
RTOT=SQRT (TOT) 
WRITE(6,105) TOT,RTOT,F1,F2,FINC 
GOTO 4 
98 FORMAT (//'FILTER LENGTH NORMALIZED FPE*// (1I7,E22.5)) 
99 FORMAT (2F9.4) 

100 FORMAT (F13.5,E13.5) 

101 FORMAT (1T2) 

102 FORMAT(/‘ENTER FILTER LENGTH (2-DIGIT INTEGER) 3", 
#* TO STOP ENTER 0;°*) 

103 FORMAT (*POHER SPECTRUM UNSCALED (INSUFFICIENT NOISE) *) 

105 FORMAT (2E13.5,3F 13.5) 

106 FORMAT (/1X,"INPUT LOWER FREQ, UPPER FREQ, FREQ', 
#"TNCREMENT;'/1X,*OUTPUTS POWER, SQRT(POWER), INPUT;* 
#/1X,'TO STOP ENTER ZEROS") 

107 FORMAT (*MEAN POWER OF THE DATA =',£10.5) 

108 FORMAT(/'ENTER 0; FOR SPECTRUM OR 1; FOR INTEGRATION") 

110 FORMAT (/* FOR SPECTRUM PLOT TYPE 1; OTHERWISE TYPE QO; '°) 

113 FORMAT (6F9.5) 

114 FORMAT("*ENTER TITLE') 

115 FORMAT (10A4) 

116 FORMAT("AUTOPICK ENTER 1; FPE LIST 23; BOTH 3; ', 
#"NEITHER 0;') 

119 FORMAT (/*FILTER LENGTH=* ,12) 

120 FORMAT (//5X,*SPECTRAL PEAKS'/*FREQUENCY POWER *, 
#* DENSITY! / (1X, F8.5,£15.5) ) 

END 
SUBROUTINE MAXENT 
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IMPLICIT COMPLEX (C) 
COMMON CX(120),CAA (120) ,CB1(120) ,CB2 (120) ,CA(120), 


#COM,CONJ, FPE (120) ,P0,PI,MM,N,ERR 


ERR=0..0 

PP=0..0 

DOTI1/I=1,N 

PP=PP+CABS (CX (I) ) **2 

PO=PP/FLOAT (N) 
FPE(1)=FLOAT (N+1) *PO/FLOAT (N-1) 
FTEMP=FPE (1) 

FPE(1) =0.0 

M=1 

CB1(1) =CX (1) 

CB2 (N-1) =CX (N) 

NN=N-1 

DO 12 I=2,NN 

CB1(1I) =CX (I) 

CB2 (I-1) =CX (I) 

COM= (0.0,0.0) 

DEN=0.0 

DO 13 I=1,NN 

CCOM=COM+CONJG (CB1(I)) *CB2 (I) 
DEN=DEN+CABS (CB1 (I) ) **2+CABS (CB2 (I) ) **2 
CA (M) =2. O*CMPLX (REAL (COM) /DEN, AIMAG (COM) /DEN) 
PO=P0* (1.0-CABS (CA (M) ) **2) 

FPE (2) =ALOG10 ( (FLOAT (N+2) *PO/FLOAT (N-2) ) /FTEMP) 
IF(PO.LE.0.0) PO=1.001 

IF (PO.LE.0.0) ERR=1.0 

DO 14 M=2,MM 

M1=M+1 

NN=M-1 

NM=N-M 

DO 15 I=1,NN 

CAA (I) =CA (I) 

COM= (0.0,0.0) 

DEN=0.0 

CONJ=CONJG (CAA (NN) ) 

DO 16 I=1,NM 

CB1(I) =CB1 (I) -CONJ*CB2 (I) 

CB2 (I) =CB2 (I+1)-CAA (NN) *CB1 (I+1) 
COM=COM+CONJG (CB1(I)) *CB2 (I) 
DEN=DEN+CABS (CB1 (I) ) ** 2#CABS (CB2 (I) ) **2 
CA (M) =2. OXCMPLX (REAL (COM) /DEN, AIMAG (COM) /DEN) 
PO=P0* (1.0-CABS (CA (M) ) **2) 

FPE (M1) =ALOG10 ( (FLOAT (N+M1) *PO/FLOAT (N-M1)) /FTEMP) 
DO) WA tenn 

CA (I) =CAA(I)-CA(M) *CONJG (CAA (M-TI) ) 
CONTINUE 

RETURN 

END 

FUNCTION POWER (F) 

COMPLEX CA (120) ,COM,CSPACE (480) ,CONJ 
COMMON CSPACE,CA,COM,CONJ,Z(120) ,PO,PI,MM 
COM=(0.0,0.0) 
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DO 19 J=1,MM 

ZZ=—-PI*F*J 

COM=COM+#CMPLX (COS (ZZ) , SIN (ZZ) ) *CA (J) 

COM=1.0-COM 

POWER=P0 /CABS (COM) **2 

RETURN 

END 

SUBROUTINE ROTATE(X,Y,ANGLE) 

SUBROUTINE ROTATE WILL ROTATE A UNIT VECTOR CLOCKWISE 
ABOUT THE X, Y, OR Z AXIS BY THE NUMBER OF RADIANS 
GIVEN IN THE INPUT PARAMETER "ANGLE". THE FIRST TWO 
PARAMETERS ARE THE VECTOR COMPONENTS WHICH WILL BE 
AFFECTED BY THE ROTATION. If IS IMPORTANT THAT THEY BE 


ENTERED IN THE CORRECT ORDER (Y FOLLOWS X, Z FOLLOWS Y, 
OR X FOLLOWS Z)~ THE MODIFIED COMPONENTS ARE STORED 
BACK INTO THESE TWO INPUT VARIABLES. 
XX=X*COS (ANGLE) —Y*SIN (ANGLE) 

Y=Y*COS (ANGLE) +X *SIN (ANGLE) 

Y=Xe 

RETURN 

END 
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APPENDIX 6 
Chi-squared Test for Log-normal Distribution 


of Power Densities 


Expected Observed (O 3-53) 
Frequency {£) Frequency (QO) 

6.45 10 a Fs SS) 
18.10 vo 0.53 
ZleDi 24 0.45 
22.67 28 Ae 25 
10.11 8 0. 44 


2244 3 0.13 


Third column total (chi-squared) = 4.75 
Degrees of freedom = 6 - 3 = 3 


Maximum "passing" value (5 percent level) = 7.82 


2/E 


132 


oo . 1% 


n Beraupe- Ei3y Istod: eo 


> Eis 3 = so beer 10 


6.0 = = thew somo 2 oniaw "page sq wea 


a) 

a r ; 7 
i 

: - oy 


~ 


: | a a 
; wey 


ahs by 


Ves ot es f 
; epee ere gly: 
r UWLar? | Rit nae 
fa ann i i | 4 
re ri : Na no ; 
1, aa 


y % : 4 ‘ 
hi . } iy : a 
wd ih 

aS at 


Searriss aete So 


Shape Se nn 


( RAGS ties Pil 
ag I Ss os Ss : ss 


Sy, 
sar hmnipsoagesine 
Aseria cates 


eeiene re ete 
SPs ~ 


Eo 


shes: ERE 
: : . ne : = - : dectgh 

qe 5 7 co : 

ae 

aay 


Dean 
LSA Sy wl AS 


ze Sees 
RANA ani ect gl en 
EU mstre Serie lee mois 


SER NRL Rit riko, 
fencpinaty ato ee ed 


lat leila 


oe. 
te 


5 ; eel IS Cystine seem oe 
: a nf B % ft . és - Dyceene nell 
peal ee ere rit Arey eens 3 Oberon ese 5 at = Wear Sey 
eds oe : ed me s ; 2 ; 5 % SE eres Sey : cs ahs ince aya se 
eatin ele eR ap - es ‘ iokrae 2 as : rh pean eee te : = SLs 
on eae, = _ - z ” ¥ < a . ; + oe pcre eta PA Nd meme 
LON hare ta ene - Z 
rnin ay tarieshiontne ep he a 


SUT aR: 


= : seer a mementos) 
oe Peat innsl 


Seep nee rane eee 
BS NR Es eee 


Dir ilvaite 
: a : Ln . : . mo ule: Cle es 
Fey enip noe Eiri eS 


